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Abstract. In this paper, we consider one dimensional adiabatic quasi-periodic Schrodinger operators 
in the regime of strong resonant tunneling. We show the emergence of a level repulsion phenomenon 
which is seen to be very naturally related to the local spectral type of the operator: the more singular 
the spectrum, the weaker the repulsion. 

■ Resume. Dans cet article, nous etudions une famille d'operateurs quasi-periodiques adiabatiques 

CD , dans un cas d'effet tunnel resonant fort. Nous voyons l'apparition d'un phenomene de repulsion de 

niveaux fort qui est relie au type spectral local de l'operateur : plus le spectre est singulier, plus la 
repulsion est faible. 

X> 

<D 

0. Introduction 

In [3], we studied the spectrum of the family of one- dimensional quasi-periodic Schrodinger 
J> ■ operators acting on L 2 (R) and defined by 



d 2 

q (0.1) Hz A = -— i/)(x) + (V(x - z) + acos{sx))ijj{x), 

q ■ where 

(HI): V : R — > R is a non constant, locally square integrable, 1-periodic function; 
. (H2): e is a small positive number chosen such that 2-k/e be irrational; 
(H3): z is a real parameter indexing the operators; 

(H4): a is a strictly positive parameter that we will keep fixed in most of the paper. 

\ To describe the energy region where we worked, consider the spectrum of the periodic Schrodinger 
g ' operator (on L 2 (R)) 

> (0.2) = + 

C3 [ We assumed that two of its spectral bands are interacting through the per- E — a E E + a 

turbation a cos i.e., that the relative position of the spectral window J-{E) := E 1 — 3 

[E — a, E + a] and the spectrum of the unperturbed operator H is that shown 
in Fig. 1. In such an energy region, the spectrum is localized near two se- 
quences of quantized energy values, say {Eq^)i and (E% (see Theorem 1.1); 
each of these sequences is "generated" by one of the ends of the neighboring pjg ure 1 ■ Bands in 
spectral bands of H . In [3], we restricted our study to neighborhoods of such i n t erac tion 
quantized energy values that were not resonant i.e. to neighborhoods of the 
points Ejp that were not "too" close to the points (Ev ^ for {/i, u} = {0, ir}. 

Already in this case, the distance between the two sequences influences the nature and location of 

the spectrum: we saw a weak level repulsion arise "due to weak resonant tunneling" . 

Similarly to what happens in the standard "double well" case (see [20, 10]), the resonant tunneling 
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begins to play an important role when the two energies, each generated by one of the quantization 
conditions, are sufficiently close to each other. 

In the present paper, we deal with the resonant case, i.e. the case when two of the "interacting 
energies" are "very" close to each other or even coincide. We find a strong relationship between the 
level repulsion and the nature of the spectrum. Recall that the latter is determined by the speed of 
decay of the solutions to (H Zj£ — E)ip = (see [9]). Expressed in this way, it is very natural that the 
two characteristics are related: the slower the decay of the generalized eigenfunctions, the larger the 
overlap between generalized eigenfunction corresponding to close energy levels, hence, the larger the 
tunneling between these levels and thus the repulsion between them. 

Let us now briefly describe the various situations we encounter. Therefore, we briefly recall the 
settings and results of [3]. Let J be an interval of energies such that, for all E G J, the spectral 
window T(E) covers the edges of two neighboring spectral bands of H and the gap located between 
them (see Fig. 1 and assumption (TIBM)). Under this assumption, consider the real and complex 
iso-energy curves associated to (0.1). Denoted respectively by and T, they are defined by 

(0.3) T R := {(«, C) G M 2 , E(«) + a ■ cos(C) = E}, 

(0.4) T := {(«, C) e C 2 , E(«) + a ■ cos(C) = E}. 

where E(k) be the dispersion relation associated to Hq (see section 1.1.2). These curves are roughly 
depicted in Fig. 2. They are periodic both in ( and k directions. 

Consider one of the periodicity cells of r^. It contains two tori. They are denoted by 70 and 7^ 
and shown in full lines. To each of them, one associates a 
phase obtained by integrating 1/2 times the fundamental 
1-form on Y along the tori; we denote the phases by $ 
and Q n respectively (see section 1.3). 
Each of the dashed lines in Fig. 2 represents a loop on 
r that connects certain connected components of T^; 
one can distinguish between the "horizontal" loops and 
the "vertical" loops. There are two special horizontal 
loops denoted by j hfi and j h y, the loop j h)0 (resp. j ht7r ) 
connects 70 to 7^ — (2ir, 0) (resp. 70 to 7^). In the 
same way, there are two special vertical loops denoted 
by Jv,o an d 7^ ; the loop j Vt0 (resp. 7„ >7r ) connects 70 
to 70 + (0, 27r) (resp. 7^ to 7^ + (0, 2ir)). To each of 
these complex loops, one associates an action obtained 
by integrating —i/2 times the fundamental 1-form on T 
along the loop. For a G {0, n} and b G {v, h}, we denote 
by S a7 b the action associated to 7^. For E real, all these actions are real 
that they all be positive. Finally, we define tunneling coefficients as 

t a ,b = e~ Sa > b/£ , a G {0,tt}, be {v,h}. 

'quantized energies" in J (see Theorem 1.1). 




Figure 2: The adiabatic phase space 
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Each of the curves 70 and j n defines a sequence of 
They satisfy the "quantization" conditions 
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denote real analytic functions of E that are small for small e. By Theorem 1.1, for 
near each Ea\ there is one exponentially small interval la such that the spectrum of 

in J is contained in the union of all these intervals. The precise description of the spectrum in 
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where o(l 

e {o,tt} 



the interval la depends on whether it intersects another such interval or not . Note that an interval 
of type 7q can only intersect intervals of type 1^ and vice versa (see section 1.3.3). 

The paper [3] was devoted to the study of the spectrum in intervals ljP that do not intersect 
any other interval. In the present paper, we consider two intervals 1^ and li l ^ that do intersect and 
describe the spectrum in the union of these two intervals. It is useful to keep in mind that this union 
is exponentially small (when e goes to 0). 

There are two main parameters controlling the spectral type: 

tvfltviT . / max(t 1 , o, t vn ) 2 
and p = 2J 

th,0 th,n y thfl th,n 

Here, as we are working inside an exponentially small interval, the point inside this interval at which 
we compute the tunneling coefficients does not really matter: over this interval, the relative variation 
of any of the actions is exponentially small. 
Noting that p > r, we distinguish three regimes: 

• t > 1, 

• p 3> 1 and r < 1, 

• p < 1. 

Here, the symbol ^> 1 (resp. <C 1) mean that the quantity is exponentially large (resp. small) in | 
as £ goes to 0, the exponential rate being arbitrary, see section 1.5. 

In each of the three cases, we consider two energies, say E and E v , satisfying respectively the first 
and the second relation in (0.5), and describe the evolution of the spectrum as Eq and E„ become 
closer to each other. As explained in Remark 1.2 in [3], this can be achieved by reducing e somewhat. 
As noted above, when moving these two energies, we can consider that the other parameters in the 
problem, mainly the tunneling coefficients, hence, the coupling constants r and p stay constant; in 
particular, we stay in one of the three cases described above when we move E and E n closer together. 

Assume we are in the case r > 1. Then, when E and E n are still "far" away from each other, 
one sees two intervals containing spectrum (sub- intervals of the corresponding intervals Iq and I n ), 
one located near each energy; they contain the same amount of spectrum i.e. the measure with 
respect to the density of states of each interval is e/2ir; and the Lyapunov exponent is positive on 
both intervals (see Fig. 5(a)). When Eq and E n approach each other, at some moment these two 
intervals merge into one, so only a single interval is seen (see Fig. 5(b)); its density of states measure 
is e/ir and the Lyapunov exponent is still positive on this interval. There is no gap separating 
the intervals of spectrum generated by the two quantization conditions . This can be interpreted 
as a consequence of the positivity of the Lyapunov exponent : the states are well localized so the 
overlapping is weak and there no level repulsion. Nevertheless, the resonance has one effect: when 
the two intervals merge, it gives rise to a sharp drop of the Lyapunov exponent (which still stays 
positive) in the middle of the interval containing spectrum. Over a distance exponentially small in 
e, the Lyapunov exponent drops by an amount of order one. 

In the other extreme, in the case p <C 1, the "starting" geometry of the spectrum is the same 
as in the previous case, namely, two well separated intervals containing each an e/2ir "part" of 
spectrum. When these intervals become sufficiently close one to another most of the spectrum on 
these intervals is absolutely continuous (even if the spectrum was singular when the intervals were 
"far" enough). As the energies E and E n approach each other, so do the intervals until they roughly 
reach an interspacing of size \fth\ during this process, the sizes of the intervals which, at the start, 

1 In the present paper, we do not discuss the gaps in the spectrum that are exponentially small with respect to the 
lengths of the two intervals Iq and 
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were roughly of order t V Q + t^ and t v n + grew to reach the order of y/th (this number is much 
larger than any of the other two as p <C 1). When Eq and E n move closer to each other, the intervals 
containing the spectrum stay "frozen" at a distance of size y/th from each other, and their sizes do 
not vary noticeably either (see Fig. 6). They start moving and changing size again when Eq and E n 
again become separated by an interspacing of size at least y/th- So, in this case, we see a very strong 
repulsion preventing the intervals of spectra from intersecting each other. The spacing between them 
is quite similar to that observed in the case of the standard double well problem (see [20, 10]). 

In the last case, when p ^> 1 and r 1, we see an intermediate behavior. For the sake of 
simplicity, let us assume that t Vj7! - 3> t v0 . Starting from the situation when E and E n are "far" 
apart, we see two intervals, one around each point E and E n and this as long as \E — E n \ ^> t V)7T 
(see Fig. 7(a)). When \Eq — E^\ becomes roughly of size t Vj7r or smaller, as in the first case, the 
Lyapunov exponent varies by an amount of order one over each of the exponentially small intervals. 
The difference is that it need not stay positive: at the edges of the two intervals that are facing 
each other, it becomes small and even can vanish. These are the edges that seem more prone to 
interaction. Now, when one moves Eq towards E n , the lacuna separating the two intervals stays open 
and starts moving with E ; it becomes roughly centered at E , stays of fixed size (of order th/t v>1T ) 
and moves along with Eq as Eq crosses E n and up to a distance roughly t v ]7r on the other side of 
E n (see Fig. 7(b)). Then, when Eq moves still further away from E v , it becomes again the center 
of some interval containing spectrum that starts moving away from the band centered at E n . We 
see that, in this case as in the case of strong repulsion, there always are two intervals separated by 
a gap; both intervals contain a e/2tt "part" of spectrum. But, now, the two intervals can become 
exponentially larger than the gap (in the case of strong interaction, the length of the gap was at 
least of the same order of magnitude as the lengths of the bands). Moreover, on both intervals, the 
Lyapunov exponent is positive near the outer edges i.e. the edges that are not facing each other; it 
can become small or even vanish on the inner edges. So, there may be some Anderson transitions 
within the intervals. We see here the effects of a weaker form of resonant tunneling and a weaker 
repulsion. 

To complete this introduction, let us note that, though in the present paper we only considered 
a perturbation given by a cosine, it is clear from our techniques that the same phenomena appear as 
long as the phase space picture is that given in Fig. 2. 

1. The results 

We now state our assumptions and results in a precise way. 

1.1. The periodic operator. This section is devoted to the description of elements of the spectral 
theory of one- dimensional periodic Schrodinger operator Hq that we need to present our results. For 
more details and proofs, we refer to section 6 and to [2, 8]. 

1.1.1. The spectrum of Hq. The spectrum of the operator H defined in (0.2) is a union of countably 
many intervals of the real axis, say [E 2n +\, E 2n+2 ] for n G N , such that 

E\ < E 2 < E 3 < E 4 . . . E 2n < E 2n +i < E 2n +2 < • • • , 
E n — > +oo, n — > +oo. 

This spectrum is purely absolutely continuous. The points (Ej)j are the eigenvalues of the self-adjoint 
operator obtained by considering the differential polynomial (0.2) acting in L 2 ([0, 2]) with periodic 
boundary conditions (see [2]). For n E N, the intervals [E 2n+1 , E 2n+2 ] are the spectral bands, and 
the intervals (E 2n , E 2n+ i), the spectral gaps. When E 2n < E 2n+ i, one says that the n-th gap is open; 

4 



when [E 2n -i, E 2n ] is separated from the rest of the spectrum by open gaps, the n-th band is said to 
be isolated. Generically all the gap are open. 

From now on, to simplify the exposition, we suppose that 
(O): all the gaps of the spectrum of Hq are open. 

1.1.2. The Bloch quasi-momentum. Let x i— ► ip(x,E) be a non trivial solution to the periodic 
Schrodinger equation H ip = Eip such that if)(x+l,E) = fiip(x,E), Vx G R, for some fi G C* 
independent of x. Such a solution is called a Bloch solution to the equation, and \i is the Floquet 
multiplier associated to ip. One may write \i = exp(ik) where, k is the Bloch quasi-momentum of the 
Bloch solution ip. 

It appears that the mapping E \— > k(E) is analytic and multi- valued; its branch points are the points 
{E n ; n G N}. They are all of "square root" type. 

The dispersion relation k t— > E(fc) is the inverse of the Bloch quasi-momentum. We refer to sec- 
tion 6.1.2 for more details on k. 

1.2. A "geometric" assumption on the energy region under study. Let us now describe the 
energy region where our study is valid. 

Recall that the spectral window T{E) is the range of the mapping (eRwE-a cos(C). 

In the sequel, J always denotes a compact interval such that, for some n G N* and for all E G J, one 
has 

(TIBM): [E 2n , E 2n+1 ] C F{E) and F{E) c}E 2n _ u E 2n+2 [. 
where T{E) is the interior of ^F(E) (see figure 1). 

Remark 1.1. As all the spectral gaps of Hq are assumed to be open, as their length tends to at 
infinity, and, as the length of the spectral bands goes to infinity at infinity, it is clear that, for any 
non vanishing a, assumption (TIBM) is satisfied in any gap at a sufficiently high energy; it suffices 
that this gap be of length smaller than 2a. 

1.3. The definitions of the phase integrals and the tunneling coefficients. We now give the 
precise definitions of the phase integrals and the tunneling coefficients introduced in the introduction. 

1.3.1. The complex momentum and its branch points. The phase integrals and the tunneling coef- 
ficients are expressed in terms of integrals of the complex momentum. Fix E in J. The complex 
momentum £ > is defined by 

(1.1) k(C) = k(E- a cos(O). 

As k, k is analytic and multi-valued. The set T defined 
in (0.4) is the graph of the function k. As the branch 
points of k are the points (Ei) i€N , the branch points of 
k satisfy 

(1.2) E - acos(C) = Ej, j G N*. 

As E is real, the set of these points is symmetric with 
respect to the real axis and to the imaginary axis, and it is Figure 3: The branch points 
27r-periodic in (. All the branch points of k lie on arccos(R). This set consists of the real axis 
and all the translates of the imaginary axis by a multiple of ir. As the branch points of the Bloch 
quasi-momentum, the branch points of k are of "square root" type. 
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Due to the symmetries, it suffices to describe the branch points in the half-strip {(; lm( > 
0, < Re( < 7r}. These branch points are described in detail in section 7.1.1 of [3]. In figure 3, 
we show some of them. The points Q being defined by (1.2), one has < ^2n < C2n+i < 7r, 
< Im C2n+2 < Im ( 2n +3 < • • • , < Im ( 2n -i < ■ ■ • < Im Ci- 

1.3.2. The contours. To define the phases and the tunneling coefficients, we introduce some integra- 
tion contours in the complex £-plane. 

These loops are shown in figures 3 and 4. The loops 70, %, 7/1,0, lh,n, %,o an d %, n are simple loops, 
respec tively, going once aroun d the intervals [-(271,(271], [bn+i, 27T-C2n+i], [~C2n+i, -C2J, [(2n,(2n+i\, 

[C2n-1) C2n-l] and [C2n+2j C2n+2] ■ 

In section 10.1 of [3], we have shown that, on each of the above loops, one can fix a continuous branch 
of the complex momentum. 

Consider T, the complex iso-energy curve defined by (0.4). Define the projection II : (£, k) G F i— > 
( G C. As on each of the loops 70, %, 7^,0, 7/1,71-, %,o and % )7T , one can fix a continuous branch of the 
complex momentum, each of these loops is the projection on the complex plane of some loop in V. 
In sections 10.6.1 and 10.6.2 of [3], we give the precise definitions of the curves 70, 7/1,0, 7/i,tt, 7^,o 
and 7„ jVr represented in figure 2 and show that they respectively project onto the curves 70, %, 7^0, 
7h,7r, lv,o and 7„ j7r . 

1.3.3. The phase integrals, the action integrals and the tunneling coefficients. The results described 
below are proved in section 10 of [3]. 



Let v G {0,7r}. To the loop j u , we associate the phase 
integral $^ defined as 



;i.s) 



$ V (E) = -<b Kd( 



where k is a branch of the complex momentum that is con- 
tinuous on 7„. The function E 1— > <!> V (E) is real analytic 
and does not vanish on J. The loop 7, is oriented so that 
<&v(E) be positive. One shows that, for all E e J, 



'M) 



& (E) < and ^(E) > 0. 



To the loop j v ,i>) we associate the vertical action integral 
S„ defined as 



(2n 7h,7r 



(2n+2 

7«,5T 
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Figure 4: The loops for the action inte- 
grals 
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S v ,u(E) 
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where k is a branch of the complex momentum that is continuous on j VjV . The function E 1— > S VjU (E) 
is real analytic and does not vanish on J. The loop 7^ is oriented so that S V>U (E) be positive. 
The vertical tunneling coefficient is defined to be 



;i.6) 



t VjV (E) = exp ( — S ViV (E) 



The index v being chosen as above, we define horizontal action integral Sh v by 



;i.7) 



S h , v (E) 



«(C) d(, 



where k is a branch of the complex momentum that is continuous on jh,v The function E i— > Sh, v {E) 
is real analytic and does not vanish on J. The loop 7^ is oriented so that Sh, v {E) be positive. 
The horizontal tunneling coefficient is defined as 

(1.8) t htV (E) = exp (--S h , v {E) 



As the cosine is even, one has 

(1.9) Sh,o{E) — Sh,n{E) and t hfl {E) = t h ^{E) . 
Finally, one defines 

(1.10) S h (E) = S h>0 (E) + Sh,Tr(E) and t h (E) = t h , (E) ■ t h , v (E). 

In (1.3), (1.5), and (1.7), only the sign of the integral depends on the choice of the branch of k; this 
sign was fixed by orienting the integration contour. 

1.4. A coarse description of the location of the spectrum in J. Henceforth, we assume that 
the assumptions (H) and (O) are satisfied and that J is a compact interval satisfying (TIBM). As 
in [3], we suppose that 

(T): 2tx ■ minmin(ImC 2 n-2(-E), lm( 2n +3{E)) > maxmax(S h (E), S V , (E), S Vy7T (E)). 

E£j E£j 

Note that (T) is verified if the spectrum of Hq has two successive bands that are sufficiently close 
to each other and sufficiently far away from the remainder of the spectrum (this can be checked 
numerically on simple examples, see section 1.8). In section 1.9 of [3], we discuss this assumption 
further. 

Define 

(1.11) S := I inf mm(S h (E), S V , (E), S v>7r (E)) > 0. 
One has 

Theorem 1.1 ([3]). Fix E* £ J. For e sufficiently small, there exists V* C C, a neighborhood of 
E*, and two real analytic functions E 1— ► <fr (E,e) and E 1— > & n (E,e), defined on V* satisfying the 
uniform asymptotics 



1.12) %(E,e) = %(E) + o(e), ^(E,e) = <f> v (E) + o(e) whene^O, 

such that, if on 
(EP(e)) v , by 



such that, if one defines two finite sequences of points in JOV*, say (Eq^)i := (Eq\e))i and (E? ' 



(1.13) -$ (E$\e) = - + 7TZ and -K(E$'\ e) = - + nl', (I, I') £ N 2 , 
el el 

then, for all real z, the spectrum of H ze in J flK is contained in the union of the intervals 

(1.14) 4° := E$ ] + [-e~ 8o/£ , e - 5 ° /e ] and if } := E^ + [-e~ So/£ , e~ &o/£ ] 
that is 



^)nJn^cfu^|j(U^) 



In the sequel, to alleviate the notations, we omit the reference to e in the functions $o and $„-. 

By (1.4) and (1.12), there exists C > such that, for e sufficiently small, the points defined in (1.13) 
satisfy 

(1.15) i e < E® - Et 1] < Ce, 

(1.16) ^e < 4° " Et 1] < Ce. 

Moreover, for v G {0, 7r}, in the interval J H the number of points E$p is of order 1/e. 

In the sequel, we refer to the points Eq (resp. E$), and, by extension, to the intervals ig^ (resp. 
I-j?) attached to them, as of type (resp. type n). 

By (1.15) and (1.16), the intervals of type (resp. n) are two by two disjoints; any interval of type 
(resp. 7r) intersects at most a single interval of type 7r (resp. 0). 

1.5. A precise description of the spectrum. As pointed out in the introduction, the present 
paper deals with the resonant case that is we consider two energies, say Ejp and E% \ that satisfy 



1.17) \EP-E^\ <2e 



2XL 



This means that the intervals I® and ij}' intersect each other. Moreover, by (1.15) and (1.16), 
these intervals stay at a distance at least C~ l e of all the other intervals of the sequences defined in 
Theorem 1.1. We now describe the spectrum of H z e in the union 1^ U l2 ■ 
To simplify the exposition, we set 

(1.18) E :=E$\ E W :=EP, 7 :=4°, and /„ := /f >. 

In the resonant case, the primary parameter controlling the location and the nature of the spectrum 
is 



(1.19) r= t vfi (E)t (E) where E= K + E 1 _ 

V th{E) 2 

As tunneling coefficients are exponentially small, one typically has either r > 1 or r C 1. We will 
give a detailed analysis of these cases. More precisely, we fix 8 T > arbitrary and assume that either 

(1.20) V-EeKHJ, S h (E) - S Vj0 (E) - S Vj7r (E) > 5 T , (case we denote by r > 1), 



or 



(1.21) ME G K H J, S h (E) - S vfi (E) - S Vj7T (E) < -S T , (case we denote by r < 1). 

The case r x 1 is more complicated (also less frequent i.e. satisfied by less energies). We discuss it 
briefly later. 

To describe our results, it is convenient to introduce the following "local variables" 

(1.22) UE ) = ^1.^-Bl where ^G{0,7r}. 



1.5.1. When r is large. Let us now assume t>1, The location of the spectrum in IoUl n is described 
by 

Theorem 1.2. Assume we are in the case of Theorem 1.1. Assume (1.20) is satisfied. Then, there 
exist Eq > and a non-negative function e i— »• f{e) tending to zero as e —> such that, for e G (0, Eq}, 
the spectrum of H Z£ in Iq U 1^ is located in two intervals Iq and I n defined by 

I = {EeI : \£ (E)\<l + f(e)} and t = {E e h : %{E)\ < 1 + /(e)}. 

If dN e (E) denotes the density of states measure of H ze , then 

(1.23) / dN £ (E)= / dN £ {E) = — if / n/, = 0, and / dN £ {E) = - otherwise. 
Jio JU 271 JiouU 71 

Moreover, the Lyapunov exponent on Iq U I n satisfies 

(1-24) Q(E,e) = S - log (ry/l + \£ (E)\ + \UE)\) + o(l), 

where o(l) — > when e — > uniformly in E, in Eq and in E n . 

By (1.23), if Iq and are disjoint, they both contain spectrum of H ze \ if not, one only knows that 
their union contains spectrum. 

Let us analyze the results of Theorem 1.2. 



2Vo(4°) -2^(4°) -2^(4°) 



(a) |£V - E'ol > max^.^o) (b) \E V - E \ < max(^, T ,t„,o) 

Figure 5: The location of the spectrum for r large 



The location of the spectrum. By (1.22), the intervals Iq and I n defined in Theorem 1.2 are respectively 
"centered" at the points Eq and E n . Their lengths are given by 

t v , (E ) ■ (1 + o(l)) and |J ff | = -r^-^—r ■ t Vt7r (E n ) ■ (1 + o(l)). 



where o(l) only depends on e and o(l) — > when e ^ 0. Depending on \E n — E \, the picture of the 
spectrum in I U I n is given by figure 5, case (a) and (b). 

The repulsion effect observed in the non resonant case (see [3], section 1.6) is negligible with respect 
to the length of the intervals Iq and 1^. 

The nature of the spectrum. In the intervals Jo and I n , according to (1.20) and (1.24), the Lyapunov 
exponent is positive. Hence, by the Ishii-Pastur-Kotani Theorem ([17]), in both Iq and the 
spectrum of H z e is singular. 
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The Lyapunov exponent Q(E,e) on the spectrum. The general formula (1.24) can be simplified in 
the following way: 

B(E,e) = Uog (ry/1 + \UE)\) + o(l) when Eel, 

and 

B(E, e) = E - log (ry/1 + \UE)\) + o(l) when £g/ . 

If l-E^ — i?o| ^ niax(t„ j7r , t^o), then the Lyapunov exponent stays essentially constant on each of the 
intervals Iq and I n . On the other hand, if \E„ — Eq\ max(t Vj7r ,t v> o), then, on these exponentially 
small intervals, the Lyapunov exponent may vary by a constant. To see this, let us take a simple 
example. Assume that t V)0 ^ tv,n, or, better said, that there exists 5 > such that 

WE e K n J, S vfi (E) > S v>7r (E) + 5. 

If Eq and i?^ coincide, then Iq C and, near the center of I n , the Lyapunov exponent assumes the 
value 

Q(E, e) = -\ogr + o(l) = - S Vj7T (E) - S vfi (E)) + o(l). 

Near the edges of /„., its value is given by 

Q(E, e) = -\ogr + ^- log(t Vj7r (E)/t Vt0 (E)) + o(l) 

7T Z7T 

= ^-(S h {E)-23 VtV {E)) + o(l). 

So the variation of the Lyapunov exponent is given by (S Vj o(E) — S Vj7r (E)) on an exponentially 
small interval. One sees a sharp drop of the Lyapunov exponent on the interval containing spectrum 
when going from the edges of I n towards E n . 

1.5.2. When r is small. We now assume that r< 1, i.e., that (1.21) is satisfied. Then, the spec- 
tral behavior depends on the value of the quantity A n (V) defined and analyzed in section 6.2, see 
Theorem 6.1. Here, we only note that A n (V) depends solely on V and on the number of the gap 
separating the two interacting bands; moreover, it can be considered as a "measure of symmetry" of 
V: taking the value 1 for "symmetric" potentials V, this quantity generically satisfies 

(1.25) A n (V) > 1. 
Below, we only consider this generic case. 

There are different possible "scenarii" for the spectral behavior when r < 1. Before describing them 
in detail, we start with a general description of the spectrum. We prove 

Theorem 1.3. Assume we are in the case of Theorem 1.1. Assume that (1.21) and (1.25) are 
satisfied. Then, there exists Eq > and a non negative function e k /(e) tending to zero when s — ► 

such that, for e g]0, Eo[, the spectrum of H Z £ in Iq U ij} ^ is contained in E(e), the set of energies E 
satisfying 

(1.26) \r%(E)UE) + 2A n (V)\ < (2 + r%(£)l + r 2 \UE)\) (1 + /(e)). 
In section 4.2, we analyze the inequality (1.26) and prove 

Proposition 1.1. For sufficiently small e, the set S(e) defined in Theorem 1.3 is the union of two 
disjoint compact intervals; both intervals are strictly contained inside the (2e _5o//£ ) -neighborhood of 
E. 
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The intervals described in Proposition 1.1 are denoted by Iq and 1^. 
We check the 

Theorem 1.4. Assume we are in the case of Proposition 1.1. IfdN £ (E) denotes the density of states 
measure of H ze , then 

e 



dN £ (E) = / dN £ (E) 
i JU 2n 

Hence, each of the intervals I® and contains some spectrum of H ZjE . This implies that, when r C 1 
and A n (V) > 1, there is a "level repulsion" or a "splitting" of resonant intervals. 

As for the nature of the spectrum, one shows the following results. The behavior of the Lyapunov 
exponent is given by 

Theorem 1.5. Assume we are in the case of Theorem 1.3. On the set E(e), the Lyapunov exponent 
of H z £ satisfies 

(1-27) e(E,e) = ±- log {t\ME)\ + + l) + o(l), 

where o(l) — > when e — * uniformly in E and in E and E n . 
For c > 0, one defines the set 

(1.28) /+ = {Ee E(e) : e\og (r^ (E)\ + \UE)\) > c} . 
Theorem 1.5 and the Ishii-Pastur-Kotani Theorem imply 

Corollary 1.1. Assume we are in the case of Theorem 1.3. For e sufficiently small, the set E(e)nl+ 
only contains singular spectrum. 

Define 

(1.29) I- = {Ee S(e) : e log (r ^ \^(E)\ + \UE)\) < -c} . 
Theorem 1.3 implies that, for sufficiently small e, the set I~ is contained in the set 

(1.30) t^(e) = {EeR: \r%(E)UE) + 2A n (V)\ <2(l + g(e))}, 

where s h- > g(e) is independent of c and satisfies the estimate g = o(l) as e —> 0. The set S ac (£) 
consists of Jo and I n , two disjoint intervals, and the distance between these intervals is greater than 
or equal to C E\Jt h {E) (see Lemma 4.12). 

Let E ac denote the absolutely continuous spectrum of H z e . One shows 

Theorem 1.6. Assume we are in the case of Theorem 1.3. Pick v e {0,7r}. There exists r/ > and 
D C (0, 1), a set of Diophantine numbers such that 

mes (D fl (0, e)) , n u s 

S: y -^J± = i + o (e~ v/£ ) when e -> 0. 

e v ' 

• for e G D sufficiently small, if I~ fl I v ^ 0, then 

mes (I v n S ac ) = mes (/„) (1 + o(l)), 

where o(l) — » w/ien £ — > uniformly in E and E n . 
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1.5.3. Possible scenarii when r is small. As in section 1.5.2, we now assume that r < 1 and A n (V) > 
1. Essentially, there are two possible cases for the location and the nature of the spectrum of H Zf£ . 
Define 



;i.3i) 



p ■-- 



max(t 1)jVr , t V) [ 



T 



max(t„ j7r , t Vj \ 



E=E 



min(^ )7r , t v>0 ) 



E=E 



Note that p > r. We only discuss the cases when p is exponentially small or exponentially large. 



1. Here, following section 5.2, we discuss the case p < 1. If \E n — E \ <C e^/th(E), E(e) is the 
union of two intervals of length roughly e-Jth ; they are separated by a gap of length roughly E\fty l 
(see figure 6); this gap is centered at the point E. The length of the intervals containing spectrum, 
as well as the length and center of the lacuna essentially do not change as the distance E w — E 
increases up to ~ eyt^E)] after that, the intervals containing spectrum begin to move away from 
each other. 

As for the nature of the spectrum, when p is exponentially small 
and \E n — Eq\ <C e^/th(E), the intervals containing spectrum 



Figure 6: When p <C 1 



are contained in the set I~; so, most of the spectrum in these 
intervals is absolutely continuous (if e satisfies the Diophantine 
condition of Theorem 1.6). 

2. Consider he case when p ^> 1. This case is analyzed in sec- 
tion 5.3. For sake of definiteness, assume that t Vj0 <C t v>n . Then, 
there exists an interval, say I n , that is asymptotically centered 
at E n and that contains spectrum. The length of this interval is of order et Vy7T {E). 

One distinguishes two cases: 

(1) if Eq belongs to I n and if the distance from Eq to the edges of /„- is of the same order of 
magnitude as the length of I n , then S(e) consists of the interval 1^ without a "gap" of length 
roughly eth(E)/t v ^(E) and containing Eq (see figure 7(b)). Moreover, the distance from Eq 
to any edge of the gap is also of order eth(E)/t v>n (E). 

(2) if Eq is outside I n and at a distance from I n at least of the same order of magnitude as the 
length of I n , then S(e) consists in the union of I n and an interval Iq (see figure 7(a)). The 
interval Iq is contained in neighborhood of Eq of size roughly e 2 th{E) /\Eq — E n \. The length 
of Iq is of size e 2 t h (E)/\E — E w \ + et v ^(E). 



e(t h /\E -E n \ +t vfl ) 



Ejr 



£t v 



Eq E-jy 



(a) When \E n - E a \ > et v ^{E) 



(b) When \E„ - E \ < et v ^(E) 



Figure 7: The locus of the spectrum when r ^ 1 and p ^> 1 



When p is exponentially large, the Lyapunov exponent may vary very quickly on the intervals con- 
taining spectrum. Consider the case \E n — E \ <C et v ^{E). For E close the gap surrounding E , 
t 2 ^q(E) is of order 1 whereas r 2 ^ 7r (£') is exponentially small. Hence, for E near the gap surrounding 
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Eq, Theorem 1.5 implies Q(E,e) = o(l). On the other hand, at the external edges of the intervals 
containing spectrum, t 2 £,o(E) is of size roughly p 2 ; this factor being exponentially large, at such 
energies, the Lyapunov exponent is positive and given by 

e(E,e) = -\ogp+o(l). 

7T 

This phenomenon is similar to that observed for r ^> 1 except that, now, the Lyapunov exponent 
sharply drops to a value that is small and that may even vanish. In fact, on most of the 
Lyapunov exponent stays positive and, the spectrum is singular (by Corollary 1.1). Moreover, near 
the lacuna surrounding Eq, neither Corollary 1.1, nor Theorem 1.6 apply. These zones are similar to 
the zones where asymptotic Anderson transitions were found in [7]. 

1.6. The model equation. The study of the spectrum of H Zy£ is reduced to the study of the finite 
difference equation (the monodromy equation, see section 2.1): 

(1.32) ^ fc+1 = M(kh + z,E)^ k , $ fc GC 2 , k e Z, 

where h = ^modl, and (z,E) — > M(z,E) is a matrix function taking values in SL(2,C) (the 
monodromy matrix, see section 2.1.1). The asymptotic of M is described in section 2.2; here we 
write down its leading term. Assume additionally that r 2 (E) > min ue ^ ^ t V)U (E) . Then, M has the 
most simplest asymptotic (see Corollary 2.1 and Remark 2.4), and, for £e/ U I n , one has 

/ r 2 g (z)g 7T (z) + 6~ l rg (z) 

(1.33) M{z,E) ~ 

y O n rg v {z) 6 n 

where 

g u = ^ + sm(27i(z - z u )), z/G {0,7r}, 
6 n is the solution to 2A n = 6 n + 6~ l in [1, +oo[, and (z v ) v( z{o^} are constants. 

The behavior of the solutions to (1.32) mimics that of those to H ZjE ip = Eijj in the sense of Theorem 
2.1 from [7]. Equation (1.32) in which the matrix is replaced with its principal term is a model 
equation of our system. All the effects we have described can be seen when analyzing this model 
equation. 

1.7. When r is of order 1. When r is of order 1, the principal term of the monodromy matrix 
asymptotics is the one described by (1.33). If r and |£ | and l^l are of order of 1, the principal 
term does not contain any asymptotic parameter. This regime is similar to that of the asymptotic 
Anderson transitions found in [7]. If at least one of the "local variables" £, V (E) becomes large, then, 
the spectrum can again be analyzed with the same precision as in sections 1.5.1 and 1.5.2. 

1.8. Numerical computations. We now turn to some numerical results showing that the multiple 
phenomena described in section 1.5 do occur. 

All these phenomena only depend on the values of the actions Sh, S Vt o, S Vj7r . We pick V to be a 
two-gap potential; for such potentials, the Bloch quasi- momentum k (see section 1.1.2) is explicitly 
given by a hyper-elliptic integral ([12, 14]). The actions then become easily computable. As the 
spectrum of H = —A + V only has two gaps, we write a (Ho) = [Ei, E 2 ] U [E 3 , £4] U [E 5 , +00 [. In 
the computations, we take the values 

E 1 = 0, E 2 = 3.8571, E 3 = 6.8571, E A = 12.1004, and E 5 = 100.7092. 

On the figure 8, we represented the part of the (a, £')-plane where the condition (TIBM) is satisfied 
for n — 1. Denote it by A. Its boundary consists of the straight lines E = E\ + a, E = E 2 + a, 
E = E3 — a and E = E4 — a. 
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The computation shows that (T) is satisfied in the whole of A. As n = 1, one has E 2n -2 = —oo. It 
suffices to check (T) for £2n+3 = (5. (T) can then be understood as a consequence of the fact that 
E5 — E4 is large. 

On figure 8, we show the zones where r and p are large and small. So, for carefully chosen a, all the 

E = E 4 - a 



E = E*-a 




E = E 2 + a 



□ 



p < 1 



p > 1, T < 1 



r > 1 



E = E 1 + a 



Figure 8: Comparing the coefficients r and p to 1 



phenomena described in section 1.5.1, 1.5.2 and 1.5.3, that is in figures 5, 6 and 7, do occur. 

1.9. The outline of the paper. The main idea of our analysis is to reduce the spectral study of (0.1) 
to the study of a finite difference equation (the monodromy equation) the coefficients of which, in 
adiabatic limit, take a simple asymptotic form. Therefore, we use the first step of a renormalization 
procedure. Such renormalization procedure, called the monodromization, was first suggested to study 
spectral properties of the finite difference equations (on the real line) with periodic coefficients in [1]. 
We have understood that this idea can be generalized to study quasi periodic systems with two 
frequencies and, in [7], we applied it to the analysis of the differential Schrodinger equation (0.1). 
In section 2, we recall the definition of the monodromy matrix and the monodromy equation for the 
quasi-periodic Schrodinger equation. Then, we get the asymptotic of the monodromy matrix in the 
adiabatic limit in the resonant case. Note that most of the technical work was already done in [3], 
where we have got this asymptotic in a general case; here, first, we analyze the error terms in the 
general formula and show that, in the resonant case, one can get for them much better estimates, 
second we show that, in the resonant case, one can simplify the principal term of the asymptotic. 
In section 3, we prove our spectral results for the case of big r. The analysis made here is quite 
standard: similar calculations can be found in [5] and [3]. 

In section 4, we prove our results for the case of small r. This is the most complicated case. Though 
that we roughly follow the analysis preformed in [7] and [3], we have to develop some new ideas to 
be able carry out rather delicate computations. This is due to two facts: first, in the case of small r, 
one observes a rich set of new spectral phenomena, and, second, one has to control simultaneously 
several objects having different orders of exponential smallness in e. 
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In section 1.5.3, we analyze the results of the previous section and describe the possible spectral 
scenarios for small r. 

The main goal of section 6 is to study A n (V), the quantity responsible for the gap between the two 
resonant intervals containing spectrum. In, particular, we show that, generically, it satisfies (1.25). 



2. The monodromy matrix 

In this section, we first recall the definitions of the monodromy matrix and of the monodromy 
equation for the quasi-periodic differential equation 

d 2 

(2.1) — —ip(x) + (V(x~z)+acos(ex))ip(x) = Eip(x), iGK, 

ax 1 

and recall how these objects related to the spectral theory of the operator H z e defined in (0.1). In 
the second part of the section, we describe a monodromy matrix for (2.1) in the resonant case. 

2.1. The monodromy matrices and the monodromy equation. We now follow [6, 7], where 
the reader can find more details, results and their proofs. 

2.1.1. The definition of the monodromy matrix. For any z fixed, let (ipj(x, z))j e { li2 } be two linearly 
independent solutions of equation (2.1). We say that they form a consistent basis if their Wronskian 
is independent of z, and, if for j £ {1,2} and all x and z, 

(2.2) ^(x,z + l) = i; j (x, z). 

As (ipj(x, z))j £ {i : 2} are solutions to equation (2.1), so are the functions ((x, z) \— > ipj(x + 2n/e,z + 
27r/e))j e { 1) 2}. Therefore, one can write 



(2.3) y (x + 2ir/E,z + 2ir/e) = M (z,E)V (x,z), V(x,z) 



i/>i(x, z) 



where M (z,E) is a 2 x 2 matrix with coefficients independent of x. The matrix M(z, E) is called the 
monodromy matrix corresponding to the basis (V'jOjeii^}- To simplify the notations, we often drop 
the E dependence when not useful. 

For any consistent basis, the monodromy matrix satisfies 

(2.4) detM(*) = l, M(z+1) = M(z), Vz. 

2.1.2. The monodromy equation and the link with the spectral theory of H z e . Set 

, \ , 2tt 

(2.5) /i=— modi. 

Let M be the monodromy matrix corresponding to the consistent basis ("0j)i=i,2- Consider the 
monodromy equation 

(2.6) F(n + 1) = M(z + nh)F(n), F(n)eC 2 , Vn e Z. 

The spectral properties of H z e defined in (0.1) are tightly related to the behavior of solutions of (2.6). 
This follows from the fact that the behavior of solutions of the monodromy equation for n — ► ±oo 
mimics the behavior of solutions of equation (0.1) for x — > =Foo, see Theorem 3.1 from [7]. 

15 



2.1.3. Relations between the equations family (0.1) and the monodromy equation. Here, we describe 
only two consequences from Theorem 3.1 from [7]; more examples will be given in the course of the 
paper. One has 

Theorem 2.1. Fix E G R. Let h be defined by (2.5). Let z \— > M(z,E) be a monodromy matrix 
for equation (2.1) corresponding to a basis of consistent solutions that are locally bounded in (x,z) 
together with their derivatives in x. 

Suppose that the monodromy equation has two linearly independent solutions (x+,X-) su °h that, for 
some C > 0, for n E N, one has ||x + (n)|| + n )\\ — Ce~ n ^ c . Then, E belongs to the resolvent 

set of the operator H z e . 

The proof of this theorem mimics the proof of Lemma 4.1 in [7]. 

The second result we now present is the relation between the Lyapunov exponents of the family of 
equations (0.1) and of the monodromy equation. 

Recall the definition of the Lyapunov exponent for a matrix cocycle. Let z \— > M{z) be an SL(C, 2)- 
valued 1-periodic function of the real variable z. Let h be a positive irrational number. The Lyapunov 
exponent for the matrix cocycle (M, h) is the limit (when it exists) 

(2.7) 0(M,h)= lim -log \\M(z + Lh) ■ M(z + (L - l)h) ■ ■ ■ M(z + h) ■ M(z)\\. 

L— >+oo L 

Actually, if M is sufficiently regular in z (say, belongs to L°°), then 9(M, h) exists for almost every 
z and does not depend on z, see e.g. [18]. 
One has 

Theorem 2.2 ([6]). Let h be defined by (2.5). Let z i— > M(z,E) be a monodromy matrix for equa- 
tion (2.1) corresponding to a basis of consistent solutions that are locally bounded in (x,z) together 
with their derivatives in x. The Lyapunov exponents Q(E,e) and 0(M(-,E),h) satisfy the relation 

(2.8) e(E,h) = ^-6(M(-,E),h). 

2.2. Asymptotics of the monodromy matrix. Consider the sequences (E^)i and (Ep)i, defined 

/>/\ 

by the quantization conditions (1.13). Let E n be one of the points [E% )i>, and let E be the point 
from the sequence (E(f')i closest to E n . Define 

^, Eft + E„ , Eft — -Err 

(2.9) E = and A 



2 2 
We assume that E and E n are resonant, i.e., that they satisfy 

(2.10) |A|<2e - ^. 

We now describe the asymptotic of the monodromy matrix for the family of equations (0.1) for 
(complex) energies E such that 

(2.11) \E-E\<Ae~^, 

We shall use the following notations: 

(1) the letter C denotes various positive constants independent of z, E, E n , E and e; 

(2) the symbol 0(fi, f 2 , ■ ■ ■ f n ) denotes functions satisfying the estimate \0(fi, f 2 , ■ ■ ■ f n )\ < 
CG/il + l/al + .-l/nl). 

(3) when writing / x g, we mean that there exists C > 1 such that C 1 \g\ < |/| < C\g\ for all 
e, C) E, E n , En in consideration. 

(4) when writing / = o(g), we mean that there exists e i— » c(e), a function such that 

16 



• \f\ < c ( £ )|#| f° r an e i Cj E n and E in consideration; 

• c(e) — > when e — > 0. 

We let 

(2.12) p(z) = e 2 " |Im ^. 
Recall that <$o is the constant defined in (1.11). One has 

Theorem 2.3. Pick E* G J. There exists V*, a neighborhood of E*, such that, for sufficiently small 
e, there exists a consistent basis of solutions of (2.1) for which the monodromy matrix (z,E) i— > 
M w (z,E) is analytic in the domain jz G C; 27r|Imz| < x K e where = {E G V*; llm^l < e}. 
Its coefficients take real values for real E and z. Fix J* C V* D R, a compact interval. If E n G J* 
and satisfies (2.10), then, in the domain 

(2.13) \z G C; 27r|Imz| < — 1 x G C; \E - E„\ < 4e"^ 



T 2 gogjr + 1 rrgo\ i fO(T 2 \g\ \g\ v p) O(rr\g\ ,p) 



one has 

(2.14) W = a ( ' yu |; ' " 

1 ; V ^ J \0 (5\g\*P, T hP ) 0(p) 

where < 5 < 5q is a constant depending only on J, a is constant in {+1, —1}, and, for v G {0,7r}, 
one has 

(2.15) g„(z, E) = Cu{E) + sin(27r(z - z v )), \g\ v = \^\ + p, 

(2.16) UE)=lu-{E-E v ). 

Furthermore, r, r, j V) 9 and z v are real constants (independent of z and E), one has 



(2.17) 



(2.18) 7i 



e-T 

(2.19) e = e n (v)(i + o(i)), 

where 6 n (V) is the positive constant depending only on n and V defined in (6.4). and 

(2.20) T h = t h {E) (1 + o(l)), f Vt0 = t v>0 (E) (1 + o(l)), %^ = t v ^(E) (1 + o(l)), 
We prove this theorem in section 2.3. 

Remark 2.1. We shall also need a more detailed description of the error terms for the coefficient 
M-Q. The proof of Theorem 2.3 (see the part of the proof in section 2.3.2 also yields 

(2.21) Ml x = o-T 2 g g„ + 0" 1 + e~ s ' £ (r 2 [0(£„6r) + O(^) + O^p) + 0(p 2 )] + 0{p)) , 
where all the error terms are analytic in E and in z. 

Remark 2.2. Being a consequence of Theorem 2.2 of [3], Theorem 2.3 stays valid if one swaps the 
indexes and it, and the quantities 9 and 1/8. One thus obtains M°, a different monodromy matrix. 
Though that all the spectral results we have announced can be obtained directly by analyzing the 
monodromy equation with the matrix M"" from Theorem 2.3, to simplify the analysis, we use from 
time to time the monodromy matrix M° instead of M n . 
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Remark 2.3. The notation r (resp. £„) in sections and 1 and the notation r (resp. £„) in the rest 
of the paper denote different quantities. They differ by a constant factor of the form 1 + o(l) (when 
e^O)! 

Let us mention a case when, preserving good error estimates, one can get a more symmetric mon- 
odromy matrix: 

Corollary 2.1. In the case of Theorem 2.3, assume that r 2 > min {T vu }. Then, there exists a 

^£{0,7r} 

consistent basis of solutions of (2.1) for which the monodromy matrix (z,E) h- > M n (z,E) is analytic 
in the domain {z G C; \lmz\ < ^} x V/ ; its coefficients take real values for real E and z. Fix 
K C V* H R compact. If E n belongs to K and satisfies (2.10), then, in the domain (2.13), one has 

(2.22) M w = a ( T2g t + B ~ 1 1°) + e~ l . (°^°\ 9 h P) ° { nW 

V 0t 9k J V °(P T \9W) 0(p) 

Remark 2.4. The principal term in (1.33) is obtained from (2.22) by replacing r, 9 and g v (more 
precisely, the quantity ^ v from the definition of g v ) by the principal terms of their asymptotics (see 
also Remark 2.3). 

Proof. Let (/i, $2) be the consistent basis for which the monodromy matrix is given by Theorem 2.3. 
Then r~ 1 / 2 fi and r 1 / 2 f 2 form a consistent basis; the corresponding monodromy matrix is 

/V-1/2 \ fr 1 ' 2 
~( rV.) ( r-V 

The asymptotic representation (2.22) follows from (2.14) and the observations that 

r~ 1 p(z) = r ^p(z) < r\g\o, and rT h p(z) = -^fa p(g) = ^^p(^) < 4t|^| w . 
Thus, Corollary 2.1 is proved. □ 

2.3. Proof of Theorem 2.3. The matrix M n is introduced in section 3.3 in [3] (M 71- was denoted by 

M u ); there, its asymptotics are described in Theorem 3.1. Under the hypothesis (2.10) and (2.11), 

we improve the estimates of the error terms in the asymptotic of the coefficient M^ 2 and simplify the 

leading terms of the asymptotics of all the coefficients of M n . This will give (2.14). 

We shall use the following notations. For (zi, • ■ ■ ,z n ) 1— > g(zx,--- ,z n ), an analytic function, we 

define 



(2.23) g*(zi, ■ ■ ■ ,z n ) = g(zi, 
For Y > 0, we let 

(2.24) T Y = e~ 2wY/£ . 



2.3.1. General asymptotic representation for the monodromy matrix M". Here, we give the asymp- 
totics of M u valid without the hypotheses (2.10) and (2.11). We check 

Proposition 2.1. Pick E* G J. There exists V*, a complex neighborhood of E*, such that, for 
sufficiently small e, the following holds. Let in (2.22), 

(2.25) Y m = ±- inf wm(S , v (E),S„,„(E)), Ym = -— sup max.(S v (E), S vn (E), S h (E)). 
2-k E£jnv t 2tt EeJnv* 
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Fix < y < Y m . There exist Y > Y M and a consistent basis of solutions of (2.1) for which the 
monodromy matrix (z,E) \—> M u (z,E) is analytic in the domain |z£C: \lmz\ < ^} x Vf. Its 
coefficients are real analytic. One has 

(2.26) M n (z, E) = P(z, E) + Q(z, E) + R(z, E), 

where 

(2>27) p ^ E) = ± (Cn(z, E)C (z, E) -S n (z, E)C (z, E)^j ^ 

, ono * A cos + 9 cos ^ cos ^ _ \ sin t^£a _ Q sin cos lo 

(2.28) Q(z,£ )= * V *»/*^ 6 ■ *o £ 

V —Osm^C^z) 9 sin -j 1 sin 

and 

(2.29) i? u (z), i* 12 (z), i? 22 (z) = 0(T h ,p(z)T Y T-\p(z)T v ), 

(2 30) - R21 ^) = COS ~f ( Th > P( Z ) T Y T h 1 iP( Z ) T v ) + COS ~ ( T h) 

+ O (pT Y , pT v T h , T%, P T V T V)7T , p T v ^TyT h ) . 

In these formulae, 

(^Z.Ol J Gt,- — - , O n — — , 

(2.32) C = ^ . 

2 

TTie functions (z,E) ^ o£q(z,E) and (z,E) i— > a n (z,E) are analytic in \z G C : |Imz| < f } x Kf; 
t/iey are 1-periodic in z and admit the asymptotics 

(2.33) a (*, 0) = 1 + T„, e 2m ^- 20 ^) + O (T Y p(z)) = 1 + (T v , p(z)) = 1 + o(l), 

0) =1 + T,, i7r [cos(2tt(z - 2^)) + f sin(2vr(z -h- z n ))} + 0(p 2 (z)7; 2 j7r , p(z)Ty) 
=l + 0(T VtV p(z)) = l + o(l). 

The functions E \— > <& V (E), E \— > T VU (E), E \— > z v (E) and E i— > T h (E), E \—> 0(E) are real analytic 
in V*; they are independent of z. In V*, one has 

(2.35) $ v {E) = $ v {E) + o{e), 

(2.36) T h {E) = t h {E){l + o(l)), T Vt „{E) = t v ,„{E) (1 + o(l)), 

where $ v and th, t v>1/ are the phase integrals and the tunneling coefficients defined in the introduction; 

(2.37) 9(E) = 9 n (V) (l + o(l)), 

where 9 n (V) is the constant defined in (6.4) (it is positive and depends only on n and V); 

(2.38) z' u (E) = 0(l/e). 

In all the above formulae, all the error terms are analytic in E and z. Finally, in the error term 
estimates, T V (E) = max{T Vj0 (E) , T V1T (E)}. 

Remark 2.5. Proposition 2.1 stays valid if one swaps the indexes 0, ir, and the quantities 9, 1/9. 
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Proposition 2.1 differs from Theorem 3.1 from [3] by more precise estimate of the coefficient R 2 i, 
and, an additional information provided by Theorem 2.2, formula (2.26) and Lemma 3.4 in [3]. So, 
to prove Proposition 2.1, we have only to check the estimate for R± 2 . 

Proof of the estimate for R 2 %. We now analyze the structure of the term R 2 \ in detail. Therefore, 
recall the description of the matrix M n given in section 3.2 of [3] (it was denoted by M u ). It has 
the form 

(2.39) = + *«W. U { .) = U { .) = \(1. J, z 



where 7(2 + h) = y ^lz) e l *' r ^ e - Here and below, we often drop the dependence on E. 

The function (z,E) 1— > a n (z,E) is analytic in G C : |Imz| < | } x V^f; it is 1-periodic in 2 and 
admits the asymptotic representations (see Theorem 2.2 and formula (2.26) in [3]) 

(2.40) a, = 1 + T v ,y^ z ~^ + O (T Y p(z)) = 1 + 0(p(z)T v , v ) = 1 + o(l). 

The branch of the square root in the definition of 7 is chosen so that j(z) = e~ l ®*/ £ (l + o(l)). 
The coefficients A n and B n are described by the formulae 



a n e^C * .i^io /0 + 1/0 T h + O* O n /6 + O* 6 O n O 



(2.41) A w = 2 + a. «S e l — 

J-h L ^ 4 L h J-h 

and 

(2.42 J B n = 2 — 1 - ot'ft e < 1 

where, for v G {0,7r}, the factor O v satisfies the estimate (see formulae (5.18) and (1.10) in [3]) 
(2.43) O v = 0{Tl T Y p(z), T h T %u p{z)). 

Note that formulae (2.41) and (2.42) are respectively the formula (5.52) and, up to a constant factor, 
the formula (5.53) in [3]. The constant factor is omitted as it can be removed by conjugation (see 
the explanations in the last lines of section 5.3.2 in [3]). 
Finally, we recall that the functions a n and a n are related by the formula 

(2-44) a^z) = a w (z) 

which is formula (3.16) in [3]. 

Now, we are ready to prove the estimate (2.30) for R 2 \. For an analytic function /, we let 3(/) = 
Representation (2.39) implies that 



2/ 



M 21 = 3(D), where D — Z ~^ ^ 



Substituting representations (2.41) and (2.42) into this formulae, we get 

M 2l = c + Ci + c 2 + c 3 + c 4 , 

where 

T h( ^ / , 2C, 



c = -6* 3 (^a e = J , c x = — is [a n a Q e * J , c 2 = — — 3 [a e * O w j , 

1„ / i(^-$ ) 1 „ / „. i(^-$ ) \ 

c 3 = -3 (^c^e * {O n /9 + O 6>)J , c 4 = —3 ^a «vre e (A-OqJ . 
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Using the representations (2.33) and (2.34) for a and a n , we get 

c = -6 sin — CV + cos — 0(pT vfi ) + 0(p T v>0 T v>7r ), 

Cl = sin T ,/4 + 0(pT w T hj p 2 ^) = cos 0(T h ) + cos ^ 0(T h ) + 0(pT v T h ), 

where, in the last step, we have also used that 0(pT v ) = o(l) for |Im£| < y < Y m . In view of the 
estimate (2.43) for O w , we get also 

c 2 = cos ^ O (T h , pT Y /T h , pT V:W ) + O (pT v , n T h , p 2 T v ^T Y /T h , p 2 T^) . 

Finally, as |oto| + \a n \ = 0(1), in view of (2.43) and as \Imz\ < y < Ym < Y, we get 

c 3 = (T 2 , P T Y , P T v T h ) , c 4 = O {Tl p 2 T 2 /T h , p 2 T 2 v T h ) . 

These estimates imply the announced representation for M^, and complete the proof of the estimate 
for R 21 . D 

2.3.2. Asymptotics of the monodromy matrix in the resonant case. We now simplify the asymptotics 
for the coefficients of the monodromy matrix given by Proposition 2.1 in the resonant case. 
In this section, J* always denotes a compact interval in V* fl J. 
For v G {0, 7r}, we let 

(2.45) a v = — sin — — — — , and a = o~o0~ w . 

Note that, by the definition of E u , one has a v G {+1, —1}, and a G {+1, — 1}. 
Define 

(2.46) f VjU = T VjU (E), f h = T h (E), and 9 = 6(E). 

Clearly, these quantities satisfy (2.20) and (2.19). 
First, for later use, we recall 

Lemma 2.1 ([3], Lemma 2.1). There exists a neighborhood of J*, say V*, and C > such that, for 
sufficiently small e, for E G V* and v G {0, ir}, one has 

\^(E)\ + \^(E)\<C, and - < \& V (E)\. 

Now, we check two simple lemmas. 

Lemma 2.2. Pick is G {0,7r}. Fix C > 0. Fore sufficiently small, if E u G J* and\E — E v \ < Ce~ 5 °^ £ , 
one has 

(2.47) cos = a„^- ■ (E - E v ) (l + O (e^e^)) 

(2.48) sin = -a„ + (^V 25 ^) . 

Proof. Both estimates follow from Lemma 2.1 and the definition of E v . We note only that to 
get (2.48) one uses (2.47). This completes the proof of Lemma 2.2. □ 
Estimate (2.47) implies that 



(2.49) cos ^ = <t v T v , v £ v (l + O (^V^)) , and 



cos 
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where £ u are the local variables defined in (2.16). 
We now prove 

Lemma 2.3. Let E e J*. For e sufficiently small, for E satisfying (2.11), one has 

T h (E) = T h (1 + 0{e~ 1 e-^)\ T vfi {E) = f v>0 (1 + O^e"^)), T v>7r (E) = f Vt1( (1 + O^V^), 

= 0(1 + 0(e-^)), z„(£) = z v {E) + O^e-^) 

Proof. Prove the representation for 7X. Recall that, in V*, a neighborhood of E* independent of 
e, Th admits the asymptotics from (2.35). This and the Cauchy estimates for the derivatives of 
analytic functions imply that, for E in any fixed compact of V*, one has logT/J < Ce" 1 . So, for 
E satisfying (2.11), we get Th(E) = IX (1 + 0(e~ 1 e~ So ' £ ). The estimates for T Vj q and T V)7T are proved 
similarly. 

Furthermore, (2.37) implies that, for E in any fixed compact of V*, one has \9'{E)\ < C. This implies 
the asymptotic representation for 9. The asymptotic representation for z u follows from (2.38) in the 
same way. This completes the proof of Lemma 2.3. □ 
We now derive simplified asymptotic representations for the functions C n , S n and Co defined in (2.31) 
and (2.32). Redefine z := z (E), and z n := z^i^E) — 2irh. We prove 

Lemma 2.4. Let E n e J*, let (2.10) hold. Fore sufficiently small, for all (z, E) in the domain (2.13), 
one has 

C = a f vfi go + e-V^T^o O (£„, p) + O (T Y p) , 

Or = <T v f v ^g v + e- l e- 5 °' e f v ^O p) + O (T Y p) , & = -a« + O {e" 5 °' e ) . 
Proof. Prove (2.50). From (2.32) and (2.33), we get 

C = cos — - h T Vt o cos I 2tt(z - z ) H — J + 0(pT Y ) 

with the "old" Zo- By means of Lemmas 2. 2, and 2.3, this yields 

Co = <7oT„,o(6 + sin(27r(z - z )) + e^e'^f^O^p) + 0(T y p) 
already with the "new" z . This result and (2.15) imply (2.50). 

The asymptotic representations for C n and S n are proved similarly; we only note that one uses the 
estimate pT V)U < e~ 5 °^ £ which follows from the definitions of p and So (as 27r|Im£| < 5q/e). This 
completes the proof of Lemma 2.4. □ 
Now, we are ready to derive Theorem 2.3 from Proposition 2.1. Begin by computing Mf x for E and 
E n e J* satisfying (2.11) and (2.10). By (2.26), we have 

(2.51) M^ = P n + Q n + Rn. 

By (2.27) and Lemmas 2.4 and 2.3, for sufficiently small e, we get 

Pn =W*^r (l + OCe-V^))- 
J-h 

(f vfi g + e-'e-^f^O (£ ,p) + O (TVp)) (T v ,„g* + e^e^^O fa,p) + O (T Y p)) . 
Now, recalling the definitions of r and g u , see (2.17) and (2.15), we get 

Pu = r 2 g g n + rW^O (£ &r,p£o,P&r,p 2 ) + O (p^- (e'^/e + (T w , w + f vfi )p + T Y p)^j 
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Note that, by (1.11), for 27r|Imz| < 5q/e, one has 

(2.52) (f h + T Ui0 + f v ^)p < Ce~ 25 °/ £ p < Ce^. 

Recall that Y > Y M (see Theorem 2.1) and note that T h > e~ 2nYM/£ , see (2.25). Let 5x = 2n(Y-Y M ). 
Then, 

(2.53) T Y /f h < Ce~ 5l/e . 
Therefore, we get 

(2.54) P n = a a n T 2 g g„ + r W^O (f of Ko, P&, p 2 ) + O (pe- {h+Sa)/e ) . 
Furthermore, by (2.28) and Lemmas 2.2 and 2.3, we have 

Q u = r 1 cos ^LZll +6 cos — cos — = O-^ao + Ofe"*/*). 

e ee 

Finally, by (2.29) and estimates (2.52) and (2.53), we get 

Rn = O (Th, p(z)T Y T-\ p(z)T v ) = 0(pe~ 5 ^ £ ,e~ So/e ). 

Fix < 5 < min(5o,5i). Substituting the estimates obtained for P n , Q n and Rn into (2.51), and 
using the notation o = <jQ<j n , we get the representation for Mfj announced in Theorem 2.3. The 
representations for the other coefficients of the matrix M 71 " are proved with the same technique. 
We only briefly comment on how to derive the representation for M^. By (2.26) and (2.27), one has 

m£ = q 21 + r 21 . 

First, one shows that 

Q21 = crerg^/r + E-'e- 6 ^ rr- 1 0(^ lP ) + 0(e-^ /e T h p). 
The representation for M 2X follows from the one for Q 2 \ and the estimate 

R 21 = e- s l e rr~ l O{pi^p 2 ) + e- s ' £ 0{T hP ) 
which follows from (2.30) and the estimates 

cos O (Th, p(z)T v ) = e^' £ rr~ l O(^), cos ^ O ( P T Y /T h ) = e~ 5 ^ e rr^O(p^), 

cos S o (T h ) = e- 1 e- s ^O(f h ), O (p 2 T v T v ^) = e-^ £ rr~ 1 0(p), 

O (p 2 ^ 1 ) = e-V'rr-^O ( P 2 ) , O (pT Y , pT v T h , T 2 ) = e^ £ 0(pT h ). 
We omit further details of the proof of Theorem 2.3. □ 

3. The case of large t 

In this section, we prove Theorem 1.2. We fix E* e J, assume that e is so small that Theorem 1.1 
holds, and systematically use notations from this theorem. 

In this section, we work under the condition r 3> 1; by (1.20), this means that, for some 5 T > 0, 
(3.1) r > e 5 ^ £ . 

Let V = {E E C : \E - E\ < 4e~ <5o/e }. Here, we study the spectrum of H Z , E in R = V n R. 

As r ^> 1, we shall use the monodromy matrix M described by Corollary 2.1. 
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3.0.3. Intervals containing spectrum. Fix N > 0. For each v G {0, n}, we let 

I V (N) = {E G J; \UE)\<l + e N }. 

We prove 

Theorem 3.1. For sufficiently small e, the spectrum of H Z£ in R is contained in Iq(N) U I n (N). 
The proof of this theorem is based upon 

Proposition 3.1 ([1], Proposition 3.1). Fix E G R and define h by (2.5). Let z i— * M(z,E) be a 
monodromy matrix for equation (2.1) corresponding to a basis of consistent solutions that are locally 
bounded in (x, z) together with their derivatives in x. 
Define 

(3.2) p(z) = M 12 {z)/M 12 {z - h), v{z) = M n (z) + p{z) M 22 {z). 
Suppose that 

(\ \ 2 

(3.3) min |Mi2(z)| > 0, max \p(z)\ < I - min \v(z)\ , indp = indv = 0, 



2 z& 

where indg is the index of a continuous periodic function g. 
Then, E is in the resolvent set of (0.1). 

Note that the proof of this proposition is based on Theorem 2.1. 

Proof of Theorem 3.1. It suffices to prove that, for e small enough, for E in R \ (I (N) U I n (N)), 
the monodromy matrix M described in Corollary 2.1 satisfy the assumptions of Proposition 3.1. 
Below, we always assume that z G R. 

1. In terms of the coefficients of M 71 " (see Corollary 2.1), define the function p by (3.2). We prove 
that, for sufficiently small e, for E G R \ Io(N) and z G R, one has 

(3.4) \M? 2 \ > 0, and \p\ < Ce~ N . 
By Corollary 2.1, we have 

(3.5) Mf 2 = arg (l + e^O(\g\ / g )). 
By (2.15), for E G R \ I (N), one has 



Mo 



9o(z ~ h) 



C £ - N (l + e N ) 



(3.6) \g \ >e N , 1 < 
For e sufficiently small, (3.5) and (3.6) imply (3.4). 

2. Here, we assume that E G R\ (Io(N) U I n (N)). In terms of the coefficients of M n , define the 
function v by (3.2). Check that 

(3.7) v = ar 2 g og7T (l + o(l)) and \v\ > Ce 2N e 2 ^ /e . 
We have 

3.8 ^oS-tt > (Co — 1)( sTtt -1) >£ and < — —— -<£ (2 + e ) . 

\9o9n\ (|?o| - lXI6r| - 1) 

Using the asymptotics of M n given by (2.22), estimates (3.8) and (3.4), we get 

v = ar 2 g o97T + 9' 1 + e"! O{r 2 \g\ \g\^ 1) + p{z){9 + 0(e"l)) = ar 2 g g n (l + O^e"!, £" 37V r- 2 )). 

In view of (3.1), this implies the first estimate in (3.7). The latter, (3.1) and the first estimate 
from (3.8) imply the second estimate in (3.7). 
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3. Steps 1 and 2 imply that, for sufficiently small e and E G R \ (Iq(N) U I n (N)), the matrix M 
satisfies the conditions of Proposition 3.1 (note that the equalities indv = indp = follow from the 
inequalities \v\, \M^ 2 \ > and the fact that the coefficients of the monodromy matrix M n are real 
valued). This implies the statement of Theorem 3.1. □ 
One has 

Corollary 3.1. Fix N > 0. For sufficiently small e, the interval Iq(N) (resp. I n (N)) is contained 
in the (Cet V fi(E)) -neighborhood (resp. (Cei Vj7T (E)) -neighborhood) of the point Eq (resp E^). 

Proof. The result follows from Theorem 3.1, the definitions of I U (N), the definitions of £„, see (2.16) 
and (2.18), Lemma 2.1 and (2.20). □ 
We shall also use a rougher result 

Corollary 3.2. For sufficiently small e, in the case of Theorem 3.1, Iq(N) (resp. In(N) ) is contained 
in the (Cse~ 2&0 I s ) -neighborhood of the point Eq (resp. E w ). 

Proof. This follows from the previous statement as, for each v G {0, 7r}, one has t VjU (E) < e -2<5 °/ £ . □ 

3.1. Computation of the integrated density of states. We now compute the increment of the 
integrated density of states on the intervals Iq and I n described in the previous subsection. We prove 

Theorem 3.2. Fix N > 0. For sufficiently small e, 

(3.9) / n £ (dE) = -, and / n £ (dE) = [ n £ (dE) = ^ if I (N) PI I n (N) = 0, 

Jl (N)Ul n (N) X Jlo(N) J'U 

where n £ (dE) denotes the density of states measure of H Z}£ . 
Proof. The proof of this theorem is based upon 

Proposition 3.2 ([3], Proposition 4.2). Pick two points a < b of the real axis. Let 7 be a continuous 
curve in C + connecting a and b. 

Assume that, for all E G 7, there is a consistent basis such that the following holds. 

• The basis solutions are locally bounded in (x,z) together with their first derivatives in x. 

• There exists V ('-/), a neighborhood of ^7 such that the monodromy matrix is continuous in 
(z, E) G K x V{i) and analytic in E G ^(7). 

• On 7, the coefficients of M , the monodromy matrix, satisfy the conditions (3.3) in which p 
and v are defined by (3.2) with h from (2.5). 

• the coefficients of M are real for real E and z. 

Then, the increment of the integrated density of states on the interval [a, b] is given by 

(3.10) / n £ (dE) = -^[ &rgv(z,E)dz 

where /| 7 denotes the increment of f when going from a to b along 7. 
The proof of Theorem 3.2 consists of the following steps. 

1. First, we assume that Io(N) fl J 7r (A r ) = and prove that J Iq / n \ n £ (dE) = e/2n. 

By (1.22), £0 is a non constant affine function of E. To apply Proposition 3.2, as 7 we choose the 

half circle 

7 = {E G C : \£o(E)\ = 1 + £ N+1 , lmE> 0}. 
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Recall that a < b denote the ends of 7. Now, [a, b] = I (N). So, I (N + 2) c]a, b[, and [a, b]r\I n (N) = 
0. So, by Theorem 3.1, a and b are in the resolvent set of H Z}£ . 
Note also that, for sufficiently small e, for E G 7, one has 

(3.11) |£ - £| < 2e~ So/£ . 

Indeed, by (2.10), one has \E — E\ < e - " 50 ^, and, for sufficiently small e, for .E G 7, we get 
|£ - £ | = eT^ol^o^)]- 1 |£ (£)| < Cst Vt0 (E) < Cee' 25o/e , 

where we have used the definitions of of 1, Lemma 2.1, (2.20) and the definitions of 5q and t Vt Q. 
By (3.11), (z,£)elx7 is in the domain (2.13), so, we can use the matrix M n from Corollary 2.1 
and its asymptotics (2.22). Define p and v in terms of M 71 by (3.2). To apply the Proposition 3.2, 
we need only to check that, along 7, condition (3.3) is satisfied. This follows from the 

Lemma 3.1. For sufficiently small e, for E G 7 and z£l, one has 

Mf 2 ^0, in dp = 0, \p\<Ce- N -\ 

v = ar 2 g g n (l + o(l)), zruit> = 0, |v| > C£ 2(Ar+1) e 2 ^ /e . 

Proof. When proving this lemma, one uses almost the same arguments as in the proof of The- 
orem 3.1 with N replaced by N + 1. The only difference is that now one deduces the equalities 
indp = indw = from the asymptotic Mf 2 = arg (l + o(l) and the asymptotics of v from (3.12). 
We omit further details. □ 
As the integrated density of states is constant outside the spectrum, formula (3.10) and the repre- 
sentation for v in (3.12) give 



2tt 2 



2 r /■! 



n £ (dE) = — J &rgv(z,E)dz 

Io(N) 



^g^gog^dz 



where, in the last step, we have used the fact that, for z and E real, the functions v, go and g n take 
real values, and, therefore, the last two integrals coinciding up to o(l), they are equal. As, for £67, 
one has \g v /£ v — 1| < 1, v G {0, 7r}, and as ^ and g v are real for real z and we get finally 



n £ (dE) = j arg(^ n )dz 

Io(N) Zn ~ Jo 



2tt 2 

7 



■ axg((E - E )(E - EJ) 



e 

27' 



where we have used the fact that only Eq is located between the ends of 7 and E n is not. 
2. Now, to complete the proof of Theorem 3.2, it suffices to check the first equality in (3.9). 
Therefore, we use the same techniques as in the previous step. So, we only outline the proof. 
Now, we take 

7 = {E G C : \E-E\ = 2e~ So/£ , lmE> 0}. 

Now, both Eq and E n are between the ends of 7. Moreover, by Corollary 3.2, both Eq and E n are 
between the ends of 7. Define p and v in terms of M 71 by (3.2). One proves 

Lemma 3.2. For sufficiently small e, for E G 7 and zel, one has 

Ml 2 ^ 0, p = 1 + o(l), iradp = 0, 1; = <tt%&{1 + o(l)), mdv = 0, |t%&t| > C/e 2 . 

Proof. One uses essentially the same analysis as when proving Theorem 3.1 and Lemma 3.1. We 
omit the details, noting only that, now, for sufficiently small e, for each v G {0, it} and E G 7, 

I > Cs-H-^E)^ - £„| > Ce-^^dE - E| - \E - E v \) > Ce^e 5 *' 5 . □ 
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This lemma immediately implies that the conditions of Proposition 3.2 are satisfied. The points £0 
and E n being between the ends of 7, one obtains 



/ n e (dE) = - ^ . aig((E - E )(E - E n )) 



5 

7T 



' Io(.N)UU(N) 

This completes the proof of (1.23). 

3.2. Computation of the Lyapunov exponent. Here, we compute the asymptotics of the Lya- 
punov exponent Q(E,e) on the intervals Io(N) and I^N) and prove 

Theorem 3.3. Fix N > 0. For sufficiently small e, for E £ Iq(N) U I n (N), one has (1.24). 

To compute Q(E, e), we use Theorem 2.2 and the matrix cocycle (M", h). In the next two subsections, 
we get an upper and a lower bound for 6(M n , h). They will coincide up to smaller order terms, and, 
thus, lead to the asymptotic formula for Q(E,e). 

In sections 3.2.1, 3.2.2 and 3.2.3, we always assume that E £ Iq(N) U I n (N). 

3.2.1. The upper bound. We now prove that 

(3.13) 6{M\ h) < 2 log (ry/l + \£ (E)\ + + C - 
Therefore, we first note that (2.22) implies that 

||Ar(z, E)\\ < CV 2 (|£ | + 1)(|^| + 1), z G R. 
Note that to get this estimate, we have used that r > 1. 

As E e (Iq(N) \Jl n (N)), then, by Theorem 3.1, for sufficiently small e, at least one of the inequalities 

|fo| < 2, and \^\ < 2. 

is satisfied. Therefore, we get 

\\M*(z,E)\\ <Cr\m + 16,1 + 1), zem. 

Now, this estimate and the definition of Lyapunov exponent for matrix cocycles (2.7) imply (3.13). 

3.2.2. The lower bound for the Lyapunov exponent. Here, we prove that 

(3.14) 6{M\ h) > 2 log (r^Jl + ME)\ + \UE)\) + 0(1). 
Therefore, we use the following construction. 

Assume that a matrix function M : C — > SL(2, C) is 1-periodic and depends on a parameter e > 0. 
One has 

Proposition 3.3. Pick e > 0. Assume that there exist y$ and y\ satisfying the inequalities < 
Ho < Hi < 00 an d such that, for any e G (0, e ) one has 

• the function z — > M(z,e) is analytic in the strip S = {z £ C; < Imz < yi/e}; 

• in the strip S = {z £ C; yo/s < Imz < yi/s} C S, M(z,e) admits the following uniform in 
S representation 

(3.15) M(z,e) = X(e)e 2mmz - ((l ^ + o(l)^ , 0, 

where X(s) and m are independent of z, and m is an integer independent of e. 
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Then, there exists E\ > such that, if < e < E\, one has 
(3.16) 6(M,h) >log|A(e)| + o(l); 

the number E\ and the error estimate in (3.16) depend only on Eq, yo. Hi and the norm of the term 
o(l) in (3.15). 

This proposition immediately follows from Proposition 10.1 from [7]. Note that the proof of the 
latter is based on the ideas of [21] generalizing Herman's argument [11]. 

Consider the case where E G Io(N). Then, one has |£o| < 2 (for e < 1). Fix < yo < y\ < where 
5 is the constant from (2.22). We shall describe the precise choice of yo and y\ later. For sufficiently 
small e and ^ < |Imz| < representation (2.22) implies that 

r + o(l)) +^e-^*—)(l + o(l)) 0^ + o(1) " 



or 



(3.17) AT = o 

We have used (3.1). Now, let y = ^ log(l + |^(£)|). Fix < A < 5/2tt. 

To compute the Lyapunov exponent for E G Iq{N) such that y < A, we choose > ^4- Then, for 
sufficiently small e, for such we get 



-27ri(2z— 2q— z-„-) 



4 



1 




oil) 



and so, Proposition 3.3 implies that 9(M, h) > log(r 2 /4) + o(l). 

On the other hand, to compute the Lyapunov exponent for all E G Iq(N) such that y > A, we choose 
yi < A. Then, for sufficiently small e, for all such E, we get 



or 
~2i 



-2iri(z — z v ) 



1 




0(1) 



and so, Proposition 3.3 implies that 9(M, h) > log(r 2 1^1/2) + o(l). 

For sufficiently small e and for E G Iq(N), the obtained lower bounds for ^(M 71 ", h) imply (3.14). 
For E G I W (N), one proves (3.14) similarly. 

3.2.3. Completing the analysis. Theorem 2.2 and estimates (3.13) and (3.14) imply that, for suffi- 
ciently small e, on I U I n , one has (1.24). This completes the proof of Theorem 3.3. 

4. The case of small t 
We now turn to the case r <C 1; by (1.21), this means that, for some 5 T > 0, 
(4.1) r < e~^ le . 

We shall assume that (1.25) holds. As before, we fix E*. G J, assume that e is so small that 
Theorem 1.1 holds, and we systematically use its notations. 
Let V = {E G C : \E - E\ < 4e" 5o/£ } and study the spectrum in R = V n 
monodromy matrix M n described by Theorem 2.3. 



Now, we use the 



4.1. The location of the spectrum. In this section, we prove Theorem 1.3. The central point of 
its proof is the a priori estimate provided by 

Lemma 4.1. Under the above conditions, for sufficiently small e, if E G RC\ a(H Z £ ), then 
(4-2) t\ME)\ + l)(\UE)\ + 1) > C(l + o(l)), 

where o(l) depends only on e, and C depends only on 6 n {V). 
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First, in sections 4.1.1 and 4.1.2, we prove Lemma 4.1. Then, in section 4.1.3, by means of this 
lemma, we get a description of the resolvent set of H Z £ inside R = {E G R satisfying (4.2)}. This 
will yield Theorem 1.3. 

4.1.1. Proof of the a priori estimate. The proof of Lemma 4.1 is based on the following construction. 
Consider the finite difference equation 

(4.3) ip(z + h) = M(z)ip(z), zER, 

where h is a fixed positive number, and M, a fixed matrix function in L°°(R, SL{2, C)). One has 
Proposition 4.1. Suppose that, for zeB, M(z) can be represented as 

(4.4) M ^ = (V °e)+M(z), 

where 8 is a real number, and this number and the matrix valued function z i— > M(z) satisfy 

(4.5) \6 + 6~ 1 \>2 and \9 - 9' 1 ] > A sup sup \Mij(z)\. 

1<M<2 zGR 

Then, there exists ip + and ip-, two vector solutions to (4.3) in L^ C (M., C 2 ), such that, for z6l, 

(4.6) det(V+(z),V-(2))>0 and + < Ce" 2 ^^, if z > 0. 

Let us first derive Lemma 4.1 from Proposition 4.1 applied to M = M 71 " and h defined by (2.5). We 
represent M n in the form (4.4) with 9 = 9. Then, for z G JR., by Theorem 2.3, one has 

|Mn| < r\\i Q \ + 1)(|^| + 1)(1 + o(l)) + 0(e"^), \M 22 \ < 0(e" 5 / £ ), 

(4.7) \M 12 M 21 \ < 9t 2 (\£ \ + + 1)(1 + o(l)) + 0(e~ s/e ). 
Only the last estimate requires to be checked. From (2.14), for z G R, we get 

\M 12 M 21 \ < (rr(\Co\ + W + o(l)) + Ce" 4 ' 6 ) (fr-Mlkl + 1)(1 + o(l)) + C^e"'/ 6 ) 

< 9t 2 (\£ \ + l)(|^l + 1)(1 + o(l)) + Ce^ e (7l,o((|£o| + 1) + ^((l&l + 1) + e"^f ft ) 
<^ 2 ((|eo| + l)(|^l + l)(l + o(l)) 

+ Ce~ 5/£ (e' 1 \E - E \+ f vfi + e' 1 \E - E n \ + f ViV + e- 5/e f h ) 

which implies (4.7). 

The estimates for the coefficients of the matrix M show that M w is similar to a matrix of the 
form (4.4) for which 

sup sup \Mij(z)\ < max{a, v^a} where a = r 2 ((|£ | + 1)(|^| + 1)(1 + o(l)) + Ce~ s/e . 
This and Proposition 4.1 imply that, if 

(4.8) 9 + 9- l >2 and a < C{9) = min {(0 - -1 )/4, (0 - ^ 1 ) 2 /(16^)} , 

then, there exist (■?/>+, "0-), the two solutions to (4.3) for M = M 71 ", that have all the properties 
described in Proposition 4.1. Define functions x± '■ Z i-^ C 2 by x±( n ) = ip±{nh + z). The functions 
X+ and X- are solutions to the monodromy equation (2.6) satisfying the conditions of Theorem 2.1. 
So, E is in the resolvent set of H z e . 

Finally, discuss the conditions (4.8). Fix < q < 1. Recall that 9 admits the asymptotics (2.37), and 
that 9 n (V) + 9-\V) = 2A n (V) > 2. Therefore, if r 2 ((|£ | + l)(|6r| + 1) < qC(8), then, for sufficiently 
small e, the conditions (4.8) are satisfied. This implies Lemma 4.1. So, to complete the proof of this 
result, we only have to check Proposition 4.1. 
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4.1.2. Proof of Proposition 4-1- Set 

(4.9) m— sup sup \Mij(z)\. 

Note that, if ip is a solution to (4.4), then e mz l h ij)(z) satisfies the same equation with M replaced 
by — M, and for a = ( 1 „], satisfies equation (4.4) with M replaced by aMa. Therefore, it 



v 1 °y 

suffices to consider the case 9 > 1. The proof then consists of five steps. 

1. We begin by the following elementary observation. Let G be a solution of the equation 

(4-10) G{z + h) = ^m±^L t ze R, 

K J y ! M 2l {z)G{z) + M 22 { Z y 

and let ip 2 (z) be a solution of the equation 

(4.11) ^ 2 (z + h) = (M 21 (z)G(z) + M 22 (z))ij 2 (z), zeR. 

Then, the vector function defined by ip(z) = ip 2 (z) i s a solution to (4.3). The proof of this 

observation being elementary, we omit it. 

2. Let 

,-»-> 1 _,/7Z^L_ 11 _ L 



2m V \ 2m 

Recall that 9 > 1 and satisfies conditions (4.5). One has 



.3-1 

,2 



(4.12) ^_ 2 ^__ 1 J, + 1=0; 

< q < 1; 

(4.13) p = ^-m(g + l)> > 1. 

The first relation is obvious; the lower bound follows from the second condition in (4.5); the upper 
bound follows from the facts that the second solution to (4.12) is greater than q and that the product 
of the solutions is equal to one; (4.13) follows from the equality 



9 _ m(9 + 1) ^ + ^(i_p_ 1 ) 2 _ 1 

and the first condition in (4.5). 

3. Let us construct a bounded solution to (4.10). For zGl and k > 1, let 

(aia\ n ( ,m M u (z)G k (z) + M 12 (z) . 

4 - 14 G k+1 {z + h)= and G {z) = 0. 

M 21 (z)G k (z) + M 22 (z) 

One has 

(4.15) sup\G k+1 (z)\ < q, keW. 

zgR 

Indeed, this estimate is valid for Go- Assume that it has been proved for a positive integer k; 
using (4.14) and (4.12), we get 

(9" 1 + m)q + m 

sup \G k+1 {z) \ < — = q. 

zgr 9 -m-mq 
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Now, check prove that {Gk(-)}^L converges in L°°. It suffices to prove that, for k G N*, 



(4.16) sup\G k+1 (z)-G k (z)\ < 



m 



p2k+l ' 



Here, p has been defined in (4.13). In view of (4.15) and (4.9), we get 

sup \Gi(z + h) — G (z + h)\ < mj (9 — m) < m/p; 

and for k > 1 and zGR 

det M(z - h) ■ (G k (z -h)- G k ^(z - h)) 



\G k+1 (z + h)-G k (z + h)\ 



(M 21 (z)G k (z) + M 22 (z))(M 21 (z)G k . 1 (z) + M 22 (z)) 
<\\G k (z-h)-G k ^(z-h)\. 

Here, we have used det M = 1; this is the only place in the proof of Proposition 4.1 where we use 
this property. The above estimates then imply (4.16). 

Denote by G the limit of {(jfc(-)}S£Lo- Clearly, G is a solution to (4.10) and satisfies |G(z)| < q. 

4. Consider equation (4.11) with G constructed in the previous step. To construct a solution to this 
equation, it suffices to define it on the interval [0, h) and continue it outside this interval by induction 
using equation (4.3); that is, for n > 1 and z G [0, h), one sets 

*fj 2 (z + nh) = (M 2l (z + (n - l)h)G{z + (n - l)h) + M 22 (z + (n - l)h))^ 2 {z + (n - l)h), 
ip 2 (z - nh) = (M 21 (z - nh)G(z - nh) + M 22 (z - nh))~ l i/j 2 (z - (n - l)h). 

Note that, for all zGl, one has \M 2l (z)G(z) + M 22 (z)\ > 9 - m(q + 1) —p > 0, and, so, for n > 1 
and z G [0, h), the second formula correctly defines ip 2 for negative z. 

Let ip 2 (z) = 1 for z G [0, h). Then, by construction, for z G [—nh, — (n — l)h) and n G N, one has 
IV^z)! < V" n ■ I n terms of ip2, we construct a vector solution to (4.3) as described in step 1. We 
denote the thus constructed solution by it satisfies 

V>_(z) = (°[ Z ^ f o r * e [0, h) and U-{-z)\\ < y/l + q ^ e - l ° sp ^ h < y/l + q 2 e^ z/h for z > 0, 

where 7 = log 0+e 2 . 

5. Construct the solution Therefore, consider the equation (4.3) with the matrix aM~ l {—z—h)a 
replacing M(z). The matrix aM~ 1 (—z — h)a can be written in the form (4.4) with the old 9 and the 
matrix aM~ 1 (—z — h)a instead of M(z). Clearly, with the matrix M, it satisfies all the hypotheses of 
Proposition 4.1. So, as when constructing we can construct ip_, a solution to the equation (4.3) 
for the matrix aM~ 1 (—z — h)a. When constructing this solution, in the last step, we normalize it 
by setting ip 2 = 1 on the interval (— h,0\. Then, we get 

$-(z) = (°[ Z ^ for z G [-h,0) and \\4>-(-z)\\ < y/l + q^ e -^ z - h)/h for z > 0. 

Here, G is a function satisfying the estimate |G(z)| < q for all z6l. 

Having constructed ip-, we define ip + : R 1— * C 2 by the formula ip+{z) = aip-(-z). The function ip + 
satisfies (4.3) for the matrix M, and one has 

i/>+{z) = (q(_ z ^ fOT z e [0, h) and U+{z)\\ < e 1 y/l + q 2 e~~< z/h for z > 0. 
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6. To complete the proof of Proposition 4.1, we need only to check that det(i/j + , 7^ 0. From 
equation (4.3), it follows that this determinant is /i-periodic. So, it suffices to consider z £ [0,h). 
Then, one has det(^+,^-) = 1 - G(z)G(-z), and so, | det(^+, V-) ~ 1| = \ G (z)G(-z)\ < q 2 < 1. 
This completes the proof of Proposition 4.1. □ 

4.1.3. Completing the proof of Theorem 1.3. The proof of Theorem 1.3 consists of three steps. In 
the first two steps, we apply Proposition 3.1 to the monodromy equations with the matrix M n and 
M° (see Remark 2.2). When we can do it, E is outside the spectrum of H Z £ . In this case, we see 
that, if E is in the spectrum of H ze , then, it satisfies (1.26). In the third step, we analyze the case 
when one can not apply Proposition 3.1; using the a priori estimate from Lemma 4.1, we see that 
then E is outside the spectrum of H Z£ . 

Below, we consider only z £ K. Fix 5^ > 0. The precise choice of this constant will be described 
later using the a priori estimate (4.2). 

1. By Theorem 2.3, one has 

M£ = <™-t(£o + sm(2n(z - z )) + e"« ■ 0(rr{\^\ + 1), 1). 
Assume that E satisfies 

(4.17) > e 5 « /e • max{l, (rT) _1 e - «}. 
Then, we get 

(4.18) M^ = arr^(l + o(l)). 
So, we have Mf 2 ^ 0. 

In terms of the monodromy matrix M 71 ", define the functions p and v by (3.2). By means of (4.18) 
and Theorem 2.3, for sufficiently small e, we get 

(4.19) p= l + o(l), 

(4.20) v = a (r 2 (£0 + sin(2vr(^ - z ))) + sin(27r(^ - z w )) + o(\^\ + 1)) + 2A n (V)) + o(l); 
when deriving the representation for u, we have also used (2.37) and (6.4). 

As the coefficients of M n are real when E and z are real, one has indp = indt? = as soon as \p\ 
is bounded away from zero, and v satisfy the second condition in (3.3). By (4.1), (4.19) and (4.20), 
there exists a function e 1— > fo(s) that is o(l) when e is small such that the first bound in (3.3) is 
satisfied when 

(4.21) \t%£ w + 2A n (V)\ > (2 + r 2 |e | + r 2 |^|) (1 + /„). 

So, for sufficiently small e, if E satisfies (4.21), it is outside the spectrum of H z e . 

2. Now, assume that E satisfies 

(4.22) \UE)\ > e s t /£ ■ max{l, rr- l e~lf h }. 

In this case, for sufficiently small e, there exists a function e 1— » /^(e) that is o(l) for 5 small such 
that, if E satisfies 

(4.23) \r%^ + 2A n (V)\ > (2 + r 2 |£o| + r 2 \^ n \) (1 + /*). 
then is outside the spectrum of -ff^, e . 

Though this result can be obtained by using directly the matrix M n , the proof becomes immediate if, 
instead of TkP, one uses the matrix M°, see Remark 2.2. Note that the conditions (4.17) and (4.22) 
are equivalent respectively to 

£o(E) > e^ E ■ max{l, 2e-^f h /f vfi } and ^(E) > e s ^ £ ■ max{l, 2e'if h /f Vjl ,} ) 
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and when one swaps the indices and n, one swaps these two conditions. So, the proof of (4.23) is 
obtained from the one of (4.21) just by swapping the indices. This completes the second step. 
3. We prove 

Lemma 4.2. Fix 8% so that 

(4.24) < < max{5 T , 6 }. 

For sufficiently small e, in the case of Theorem 1.3, all the energies E that satisfy neither (4.17) 
nor (4.22) are outside the spectrum of H z e . 

Proof. Pick E that satisfies neither (4.17) nor (4.22). Then, one has 

r 2 (\UE)\ + l)(|£o(£)| + 1) < Cr 2 e 25 ^ • (l + rr^e^f h ) (l + (rr)" 1 ^) 

= Ce^'e . (> + e -f (f^ + f Bi0 ) + e^T h ) < Ce 2 ^ ■ (e~* + e"^) , 

where, in the second step, we have used (2.17), and, in the last step, we have used (2.20) and (1.11). 
In view of the last computation and (4.24), we get that r 2 (|^ 7r (£')| + l)(|£ (-£/)| + 1) = o(l), and, in 
view of Lemma 4.1, for sufficiently small e, this implies, that E is outside the spectrum of H z e . This 
completes the proof of Lemma 4.2. □ 
Now, let f(e) = max{f (e), f^e)}. Clearly, for sufficiently small e, the function / is well defined 
and satisfies / = o(l) near 0. By Lemma 4.2 and the first two steps, we see that, if E in R is in the 
spectrum of H Z)£ , then, it satisfies |t 2 £o£tt + 2A n (V)| < (2 + r 2 |£ | + i" 2 |£tt|) (1 + /)• This completes 
the proof of Theorem 1.3. □ 

4.2. Properties of the set defined by (1.26). We now analyze (1.26) and prove Proposition 1.1. 
Let 

(4.25) G(E) = r%(E)UE) + 2A n (V), and F(E) = (2 + r%(£)| + r 2 \UE)\) (1 + /(e)), 
where the function e ^ f(e) is defined in Theorem 1.3. So, the condition (1.26) takes the form: 

(4.26) \G(E)\<F(E). 

4.2.1. Individual properties of F and G. We first discuss various properties of the functions G and 
F without comparing their graphs, i.e., without analyzing (4.26). 

Lemma 4.3. One has 

(1) G is a quadratic polynomial in E; 

(2) it reaches its maximum at E; 

(3) G(E ) = G(E 7r ) = 2A n {V) > 2; 

(4) F is affine except at the points E and E^; 

(5) F is strictly increasing for E > max{i?o, E n }, strictly decreasing for E < mm{Eo, E^}; 

(6) between E and E^, the absolute value of the derivative of F is smaller than it is for E > 
ma.x{E , E-x} or E < min{£'o, E n }; 

(7) F(^) = (2 + r 2 |e (^)|)(l + /) andF(£ ) = (2 + T 2 |U£o)|)(l + /). 



Proof. Lemma follow from (2.16) and (4.25). To prove point (2) one also uses (2.18) and (1.4). □ 
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4.2.2. The intervals I in and I ont . Now, we begin the analysis of condition (4.26). First, we describe 
the set {E E R; G(E) > F(E)}. 

Corollary 4.1. For sufficiently small e, the set I[ n = {E E R; G(E) > F(E)} is a compact interval 
of positive length. It is located strictly between the zeros of the polynomial G. 

Proof. Lemma 4.3 implies that G is concave (points (1) and (2)). As F is positive, J in (if not empty) 
is located between the zeros of G. By the points (4) - (6) of Lemma 4.3, F is convex. So, now, it 
suffices to prove, that for sufficiently small e, there exists a point E, where G(E) > F(E). Therefore, 
note that = |£o(-Eo)| < l&rCEb)!, and that = ^(E^)] < \£ (E n )\. This implies that, between E 
and E n , there is E, a point where £o(E) = £ n (E). Denote this common value by £. At E, one has 

G(E) - F(E) = r 2 |£| 2 + 2A n (V) - 2(1 + r 2 |£|)(l + /) 

= r 2 (KI - 1 - f? + 2{hjy) ~ 1 " /) " (1 + f? 
>2(A n (V)-l-/)-r 2 (l + /) 2 . 

Therefore, for sufficiently small e, one has G(E) — F(E) > A n (V) — 1 > 0. This completes the proof 
of Corollary 4.1. □ 
We now prove 

Corollary 4.2. The set I ont = {E e R; G(E) > —F(E)} is a compact interval of positive length. 
Moreover, J in is contained in I ont , the interior of I out . 

Proof. Let e be so small that I- m exists and that F > for all E G R. On the compact interval 
bounded by E n and E Q , the function G is positive (see points (2) and (3) of Lemma 4.3), whereas 
— F is negative. This, the facts that G is a concave quadratic polynomial and — F is piecewise affine 
and concave (by Lemma 4.3) imply that J out is a compact interval of positive length. The inclusion 
4i C /out follows from the inequality G{E) > F(E) > -F(E) valid for E e I m . □ 

4.2.3. The intervals I\ and I r . Let be the set where condition (4.26) is satisfied. One has 
£(e) = /out \ /in- Corollaries 4.1 and 4.2 imply 

Corollary 4.3. For sufficiently small e, the set E(e) consists of two disjoint compact intervals of 
positive length. 

We denote these intervals by I\ and / r so that I\ be to the left of / r . We finally check 

Lemma 4.4. For e sufficiently small, both l[ and I r are inside the (2e~ 5 °/ £ ) -neighborhood of E. 

Proof. It suffices to check that outside the (2e -<5 °/ £ )-neighborhood of E, one has > F(E). 

By (1.11) and Lemma 2.1, for E such that \E — E\ > 2e~ s °^ e and v E {0, 7r}, one has 



\&„(E)(E - E v )\ e-f 2 A\%{E)%{E){E - E Q ){E - E^)\ 1 
\UE)\ = ^- r|&(S)US)l = >^ 

For sufficiently small e, these estimates imply that, for \E — E\ > 2e~ s °/ £ and v E {0, 7r}, 

|G(E)| - F(E) > t*(\UE)\ - 1 - f)(\UE)\ - 1 - /) - 2(A n (T/) - 1 - /) - r 2 (l + ff > 



Ce 2 ' 

This completes the proof of Lemma 4.4. □ 



Corollary 4.3 and Lemma 4.4 prove Proposition 1.1. 
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4.3. Computing the density of states. We now compute the increments of the integrated density 
of states on each of the intervals defined in Proposition 1.1 and, thus, prove Theorem 1.4. 
We assume that 



(4.27) 




The complementary case is treated similarly, but instead of working with M 71 ", one uses the mon- 
odromy matrix M° mentioned in Remark 2.2. For sake of definiteness, we assume that o = 1 in (2.14) 
and that 



The cases where E n > E or a = — 1 are analyzed similarly. 

Our main tool is Proposition 3.2. First, we compute the increment of the integrated density of states 
on the whole set S e and, then, we compute it on one of the intervals of this set. 

4.3.1. The increment of the IDS on the set E £ . As the curve used to apply Theorem 3.2, we choose 



Lemma 4.4 implies that the set E e is strictly between the ends of 7. 

The analysis of the increment of the integrated density of states between the ends of 7 is standard, 
see section 3.1. We omit details and note only that, first, in terms of M n , one defines v and p by (3.2), 
and, then, one checks that for 7 defined in (4.29), the statements of Lemma 3.2 hold. As in part 2 
of section 3.1, this implies that the increment of the integrated density of states between the ends of 
7 is equal to e/ir. 

4.3.2. The increment of the IDS on one of the intervals o/S e . To complete the proof of Theorem 1.4, 
we pick E*, a point in between the connected the components of the set S e , and prove that, between 
E* and a point located on R outside the (2e _<5o//£ )-neighborhood of E, the increment of the integrated 
density of is equal to ej (2tt). 
We define E* = E — T h — T v>0 and 

7_i_ = {E G C; KeE = E*, E is between R and 7}. 

Now, let 7* be the curve going from E* along 7j_ to 7 and, then, along 7 to K in the clockwise 
direction. This is the curve we use to apply Theorem 3.2. 
Let us study v and p on the curve 7j_. We shall use 

Lemma 4.5. For sufficiently small e, for z G R and E = E*, one has 



(4.28) 



E n < Eq. 



(4.29) 



7 = {E G C; \E-E\ = 2 e~ 5o/e , lmE> 0}. 




Lemma 4.6. For sufficiently small e, for z G K and E G 7_i_ ; one has 
(4.33) Aff 2 ^ 0, p=l + o(l), mdp = 0. 
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Proof. By (2.14), for z G R, we have Mf 2 = rr{ (l + o(l)) +0(rr) +o(l). Note that estimates (4.31) 
are valid as at E = E* so on the whole curve r )±. Therefore, under the conditions of Lemma 4.6, 
M^ 2 = rr£o(l + o(l)). This implies all the statements of Lemma 4.6. □ 

Now, we turn to the function v. We prove 

Proposition 4.2. For sufficiently small e, for z G R, one has 

• if E G 7_l, then Rev > 2A n + o(l); 

• ifE = E„, thenv = (r 2 £ ^ + 2A n )(l + o(l)). 

Remark 4.1. When proving Proposition 4.2, we shall see that, for sufficiently small e and = £7*, 

t- 2 ^o = o(1), r 2 ^ = o(r 2 ^^ ), and r 2 ^ > o(l), 

see (4.36) and (4.39). This implies that, for sufficiently small e, the point E* is inside the interval 
L m , i.e., between the intervals of the set E e described by (1.26). 

Proposition 4.2 immediately follows from the next two lemmas. 

Lemma 4.7. For sufficiently small e, for z G K and i? e 71, one has 

(4.34) Rev = r 2 Re£ Re«^(l + o(l)) - r 2 Im£ Im^(l + o(l)) + 2A n + o(l). 

Proof. Recall that r is small. Using (2.21), the representation MJ 2 = 8 + o(l) following from (2.14), 
and the asymptotics p = 1 + o(l), see Lemma 4.5, we get for z6K and E G 7_i_ 

Re v = r 2 (Re £ Re ^ - Im £ Im + 2A n 

(4.35) + r 2 Re (f ) sin(27r(2 - z )) + r 2 Re sin(2vr(2; - ^)) 

+ e- s ' e t 2 Re (O(^) + O(fo) + O(^)) + o(l), 

where the terms 0(^0^), O(£o) anc ^ ^(^) are analytic in E 1 . Let us study the terms in the second 
and the third lines of this formula. 
First, we prove that 

(4.36) r 2 Re(£ ) = o(l), r 2 Re (&) = o(r 2 Re^ Re£ ). 
Clearly, 

Re (£ ) =£o(£*), Re (&) 
As |Co(-E*)| > 1 and r 2 |£ (£*)| < 1, see (4.31) and (4.32), this implies (4.36). 
Now, prove the estimate 

(4.37) ReO(£o^) = O(Re£ Re^) + O(Im£ Im^). 

Note that as in (2.21), the term O(£o£?r) is a real analytic function of E bounded by C|£o£tt| uniformly 
111 the (4e « ) -neighborhood of E. Therefore, it can be represented in the form 

where g is a real analytic function of E satisfying the estimates 

\g(E)\ < C for \E - E\ < 4e"^; 
\g\E)\ < Ce^, for \E - E\ < 3e~^; 
\lmg(E)\ < C|Im£|e^, for \E - E\ < 3e"^. 

As 

\BeO(£o£*)\ < \Re(^)Reg\ + |Im fo£ v )lmg\, 
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the estimates for Img and \g\ imply that 

ReO(^n) = 0(Re«£o Re^) + O(Im£ Im^) + e So/e (O(Imf Re^ JmE) + 0(lm£ n Re«£ JmE)). 
Using the definitions of £ an d £tt, on j±, we get 

Im£ Re^ Im E = Im^ im ^7r(-E'* — E n ) and Im^ Re£o JmE = lm^ lm^ 7r (E^ — E Q ). 

As both \E* — E \ and \E* — E w \ are bounded by Ce~ 5 °^ £ , we finally get (4.37). 
Using the same techniques, for each v G {0, 7r}, one also proves that 

ReO(^) = 0(Ref„) + e So/e O(lm^lmE) 

= 0(Re&) + T t , iM e^o(Im£ Im^) = 0(Re£„) + o(Im^Im^), 

where \i is the index complementary to v in {0, tt}. In view of (4.36), this implies that 

(4.38) r 2 Re(O(£ ) + O(^)) = o(l) + o(r 2 Ref Re&r) + o(r 2 lm^ Im^). 

Substituting estimates (4.36) - (4.38) into (4.35), we come to (4.34). This completes the proof of 
Lemma 4.7. □ 

Lemma 4.8. For sufficiently small e and E G 7_|_ ; one has 

(4.39) Im£ Im^<0, r 2 Re^ Re^ > o(l). 
Proof. It follows from the definitions of £o and ^ that 

(4.40) l^ Q l^= ^^ E) (Im£) 2 , r^ Re^ = A ^ E ^ (g . - E )(E. - E n ). 

For each v G {0,7r}, one has &' V {E) = &' U (E) + o(l) (which follows from (1.12) and the Cauchy 
estimates). So, in view of (1.4), we get <b' (E)<b' n (E) < 0. Therefore, (4.40) implies that Im£ Im^ < 
0, and that, if E n < E* < Eq, r 2 Re^ Re^7r > 0. To complete the proof of Lemma 4.8, we need only 
to check that, if E* < E n < E , one has r 2 Re^oRe^7r = o(l). But, for such values of E n , the second 
formula in (4.40) and the definition of E* imply that 

%(TlJT h + 2T vfi + T h ) < % 



|Re£ Re^| < ~ (T V)0 + T h f = - (f 2 /f h + 2T vfi + T h ) < — (r 2 + 2T V>0 + T h ) 

£ J-h £ 



which implies the needed estimate. This completes the proof of Lemma 4.8. □ 
Now, we are ready to prove 

Lemma 4.9. For sufficiently small e, the matrix M 71 " and the curve 7* satisfy the assumptions of 
Theorem 3.2. One has 



(4.41) / argu(x, E)dx 

Jo 



= avgG(E) 

7* 



where G(E) = r 2 £ &r + 2A n . 

7* 

Proof. The assumptions of Proposition 3.2 are satisfied as 

• as we have already mentioned, on 7* n 7, the statement of Lemma 3.2 holds; 

• on 7* fl 7_i_, the function p is described by Lemma 4.6 and, in view of the first point of 
Proposition 4.2, one has mdv = 0, and \v\/2 > A n + o(l) (recall that A n > 1 is a constant 
depending only on V and n, and that we consider the case where A n > 1). 

Now, let us prove (4.41). As both v and G are real analytic, the left and the right hand sides 
of (4.41) coincide modulo n. So, it suffices to prove this equality up to o(l). This follows from the 
observations: 
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in view of Lemma 3.2, on 7* fl 7, one has v(x, E) = G(E){1 + o(l)); 

in view the previous point and the second point of Proposition 4.2, at the ends of 7*fl7_i_ = 7^, 
one also has v(x, E) = G(E)(1 + o(l)); 

in view of the first point of Proposition 4.2 and as v is real analytic, one has 



— 7r/2 < / &rgv(x, E)dx 







< tt/2. 



7_ 



• in view of Lemma 4.8, on 7j_, KeG(E) > 2A n > 0, and so, as G is real analytic, one has 

-tt/2 < argG(£)| 7± < tt/2. 

This completes the proof of Lemma 4.9. □ 

Now, we note that argG(E)\ = —n. This follows from Remark 4.1 and the statement of Lem- 
ma 4.2.3 saying that I\ n is located strictly between the zeros of the polynomial G. So, Theorem 1.4 
follows from (3.10) and (4.41). □ 

4.4. Computing the Lyapunov exponent. We now prove Theorem 1.5. The computations are 
essentially the same as in section 3.2, but, instead of working with the matrix cocycle (M n ,h), we 
pass to an auxiliary one. 

Below, we always assume that E G S s , i.e., that it satisfies (1.26). Note that this implies that 

(4.42) T 2 ma*\ < c(t 2 \u +^1 + 1), 

Fix a constant S satisfying 

(4.43) 0< 5s <mm(5,5 T ,5 ). 

Then, by Lemma 4.2, for any E satisfying (1.26), one has either (4.17) or (4.22). We prove The- 
orem 1.5 assuming that E satisfies (4.17). The other case is treated similarly, but, instead of the 
matrix M n , one uses the monodromy matrix M° introduced in Remark 2.2. 

4.4.1. Auxiliary matrix cocycle. We use the following 

Lemma 4.10 ([3], Lemma 4.7). Let M G L°°(R, SL(2, C)) be 1-periodic, and let h be an irrational 
number. Assume that there exists A > 1 such that 

(AAA) Vx G R, A' 1 < M l2 (x) < A. 

In terms of M and h, construct v and p by formulae (3.2). Set 



(A 45) N(x) = ( ^ x )/v^) -\fP@, 

Then, the Lyapunov exponents for the matrix cocycles (M, h) and (N, h) are related by the formula 

(4.46) 9(M,h) = 6(N,h). 

Now, for M = M n and h defined by (2.5), we construct iV by formula (4.45). Under the condi- 
tion (4.17), Theorem 2.3 implies that 

(4.47) Mf 2 = arr^il + 0(pe~^ /£ )). 
So, we are in the case of Lemma 4.10. 
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4.4.2. The coefficients of N. Here, we check 

Lemma 4.11. For sufficiently small e, if E satisfies (4.17) and if 27r\lmz\ < one has 

(4.48) N= (^(t sin(2n(z -,.))+ ^ sin(2^z - z ))) 0^ + (r 2 (|£ | + |^|) + l) . 

Proof. For sufficiently small e, for z and E as in Lemma 4.11, representation (4.47) implies that 
p — 1 + o(l). This, (2.14), (4.43) and (4.42) imply that (for sufficiently small e, for z and E we 
consider) 

u = ar 2 (Z sin(2vr( 2 - ^)) + £ ff sin(27r(z - z ))) + O(r 2 (|£ | + |^l) + 1)- 
This representation and the representation p = 1 + o(l) imply (4.48). □ 

4.4.3. ^4n upper bound for Q(N, h). Lemma 4.11 implies that, for sufficiently e, for z6l and for E 
satisfying (4.17), one has 

\\N\\ <C(r 2 (|£ | + |^|) + l). 
This and the definition of the Lyapunov exponent for a matrix cocycle imply that 

(4.49) G(N, h) < log(r 2 (|£ | + U) + 1) + C. 

AAA. A lower bound for Q(N,h). Let us now show that, for sufficiently small e, for all E satisfy- 
ing (4.17) and such that 

(4-50) T 2 (|£ | + |&r|) > 1, 

one has 

(4.51) G(A^, ^) > log(r 2 (|e | + + c. 
The proof consists of two steps. 

1. For sufficiently small e, for all E satisfying (4.17) and (4.50), one has either 

(4.52) |Co| > C/r 2 and |^| < C, 
or 

(4.53) > C/t 2 and |f | < C. 

Indeed, (4.42) and (4.50) imply that |£o||£tt| < C(lCo| + 1^1), and, for sufficiently small r, this 
inequality and (4.50) imply the above alternative. 

Below, we consider only the case of (4.52). The second case is treated similarly. 

2. Fix y so that < y < 5^. Lemma 4.11 implies that, for sufficiently e, for y < 27re|Imz| < 5^ and 
for all E satisfying (4.17) and (4.52), one has 



(TT 2 

(4.54) N = —^e- 2ni{z - z ^ 

This, Lemma 3.3 and (4.52) imply (4.51). 



1 




(1) 



4.4.5. Completing the proof. From (4.49) and (4.51), we conclude that 6(N, h) = log(r 2 (|£o| + l&rl) + 
1) + 0(1). This and Theorem 2.2 imply (1.27). This completes the proof of Theorem 1.5. □ 

4.5. Absolutely continuous spectrum. Here, we prove Theorem 1.6. The proof consists of two 
main steps. As when computing the density of states, we work under the assumptions (4.27) 
and (4.28). 
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4.5.1. Properties of the set I c . Let us discuss properties of the sets I c and £ ac (e) defined in (1.29) 
and (1.30). Recall that A = \E n - E \/2. Let 

e 2 f h 



(4.55) M(A) = min |eVT h , 
First, we check 

Lemma 4.12. The set S ac (e) consists of two disjoint intervals Iq and I n such that Iq is to the right 
of Eq, I n is to the left of E n ; for each v G {0,7r} ; one has 

(4.56) \I V \ x M(A) and dist (/„, E v ) x M(A), 
where \I\ denotes the length of an interval I. 

Proof. Using definitions of £ an d see (2-16), and the ones of £ ac , see (1.30), one obtains 
Sac(^) = lo U h where 

I 7T =[E-b,E-a], I =[E + a,E + b), 



A^ + e*T hl (A n -1-g), b = y A 2 + e 2 T hl (A n + l + g), 

Note that, by (1.4) and (1.12), one obtains 7x1; furthermore, recall that A n > 1. 
The above formulae already imply that I and 1^ are disjoint and their length satisfy (4.56). Moreover, 
as Eq = E + A and E n = E — A, see (4.28), they also imply the statements on the positions of the 
intervals Iq and This completes the proof of Lemma 4.12. □ 

We shall use 

Corollary 4.4. Let d = ^ dist (E , I ). On the do-neighborhood of Iq, 
Let d % — | dist (E n , 1^). On the d n -neighborhood of I n , 

Proof. Corollary 4.4 follows from Lemma 4.12 and the definitions of £o and £„.. □ 
The following property of I~ plays important role: 

Lemma 4.13. Fix c > and v G {0,7r}. There is a positive constant C such that, for sufficiently 
small e the following holds. If I~ PI I v ^ 0, then, for all E G I u , one has 

(4.57) r 2 (|6(^)| + |^)|)<Ce- c / £ . 

Proof. Consider the case when I~ D Iq ^ 0- The complementary case is analyzed similarly. 
The statement of Lemma follows from the observations: 

• at a point of J , one has t 2 (\£ (E)\ + \£*(E)\) < e~ c/£ (as I~ n J + 0); 

• simultaneously for all E G Jo, one has |^ (-S)| + |^7r(-E , )| x Q where q is positive and independent 
of E (Corollary 4.4 gives q = Mm + A +^ A) ). 

The proof of Lemma 4.13 is completed. □ 
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4.5.2. The monodromy matrix for E e I~ . When proving our results on the absolutely continuous 
spectrum, we use Ishii-Pastur-Kotani Theorem, i.e., we control the Lyapunov exponent for the family 
of equations (0.1) using the matrix cocycle N defined in Lemma 4.10. We now study the matrix N 
constructed by formula (4.45) in terms of the matrix M n . 
Fix < 5 ac < min{c, 5, 5 T }. We prove 

Proposition 4.3. Fix y > 0. Pick v e {0,7r}. For sufficiently small e, the following holds. If 
I~ H I v 7^ 0, then, for E in the d u -neighborhood of I u and for z in the strip \lmz\ < y, one has 

(4.58) N = Y" + l/dn) ~o) + 0{e~ 5 ^). 

Proof. The cases where v — and v = ir are treated in one and the same way. We only consider the 
case of v = 0. Below, we always assume that E and described in Proposition 4.3. Construct 

the functions p and v in terms of the matrix M n . It suffices to prove that 

(4.59) p = l + 0(e- s - /e ), v = a(T%^ + 9 n +l/6 n ) + 0(e~ s ^ £ ). 
The asymptotic representation for p follows from (2.14) and the estimates 

(4.60) |fo| > Ce 5ac/£ and \rr£ \ > 1. 
Let us prove these estimates. By Corollary 4.4, we get 

(4.61) M>C»®. and rr M > M 

We consider two cases. First, we assume that A < e\fTh. Then, M(A) = e^/Th, and we get the 
estimates (using also (4.27)), 

Ifol > > Cr- 1 > Ce 5 ^ £ and rr|Co| > 



which imply (4.60). Now, assume that A > ey/T^. Then, M(A) = e 2 Tf l /A, and we get the estimates 

fT Cf 
(4-62) ^- C T~A ^ rr ^°l - X' 

In view of (2.10), the second of these estimates implies the second of the estimates in (4.60). Now, 
note that Lemma 4.13 and Corollary 4.4 imply that 



e c/£ < CCr 2 !^!)- 1 < C 



eT h 
T.nA' 



Therefore, the first of the estimates in (4.62) implies the first estimate in (4.60). This completes the 
proof of the asymptotic representation for p. 

Prove the asymptotics for v. The representation (2.14) and the already proved representation for p 
imply that 

(4.63) v = ar%Cn + n + l/9 n + O {e~ 5 ^(r 2 \^ w \ + 1) + r 2 (|£ | + |^D) ■ 

Now, note that, for all E G Io, one has (4.57) and t 2 |£o£tt| < C (as I C S ac ). Therefore, (4.63) 
imply the representation for v from (4.59). This completes the proof of Proposition 4.3. □ 
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4.5.3. Completing the proof. Using the representation (4.58), one completes the proof of Theorem 1.6 
as in section 4.5 of [3]: using standard KAM techniques, see [7], section 11, one proves that, under the 
conditions of Theorem 1.6 the Lyapunov exponent for the matrix cocycle (N, h) vanishes on the most 
of the interval I u , i.e., the interval in the neighborhood of which the representation (4.58) holds. By 
Lemma 4.10 and Theorem 2.2, this implies that the Lyapunov for the equation family (0.1) vanishes 
on most of this interval. Then, by the Ishii-Pastur-Kotani Theorem, this implies that most of this 
interval is covered by absolutely continuous spectrum. 

To complete the proof, we first transform the problem to a form suited to apply the standard KAM 
techniques, then, we complete the analysis as in [3]. The last part being standard, we only outline 
it. 

Step 1: New parametrization. One introduces <p \— > E(ip), the multivalued analytic function defined 
by the relation 

2 cosy? = Tr( leading term of N) = o-T 2 £ {E)£ n {E) + 6 n + l/6 n . 

One picks four positive constants q\, q2, c\ < 1 and C2 < 1 — c\ and considers the function (p i— > E(<p) 
on V(ci, C2), the complex (c2)-neighborhood of the interval [—1 +c±, 1 — c\\. Using (2.16), one proves 
that, if Ci and C2 are sufficiently small, then, for sufficiently small e, 

• there exists a real analytic branch of the function ip 1— > E(ip) that maps the interval [— 1 + 
ci, 1 — ci] into I v so that \I U \ — 1 + c\, 1 — Ci))| < giM(A); 

• it conformally maps V(ci,C2) into the g 2 ^(A)-neighborhood of l v \ 

From now on, one considers A^ as a function of (z, if) G {z G C : |Im z\ < y} x V(ci, C2). 
Step 2. Diagonalization of the leading term of N . Let 



(4.64) A> = 1 S" 1 A^ 1 S where S 



Clearly, together with N, the function A" is unimodular, 1-periodic in z and analytic in (z, ip) G 
{z G C : |Imz| < y} x U(ci,C2) (by the previous step, provided that c\ and C2 are chosen small 

enough). Furthermore, as A^ is real analytic, iV has the form ( ^ ^ ] (see (2.23)). Finally, for 



(z, <f) G {z G C : \lmz\ < y} x "K(ci,c 2 ), one has 



(4.65) N = ^ e °^J + 0(A) where A = e~ s ^ £ 

(by Proposition 4.3 and the previous step, provided that c\ and C2 are chosen small enough). 
3. Standard KAM theory result. Clearly, 8(N,h) = 8(N,h). Applying a standard KAM theory 
construction from section 4.5.2 of [3] to the cocycle (A^, h), one proves that, if the number h satisfies 
a Diophantine condition (see Proposition 4.4 in [3]), and A is sufficiently small, then #(A^, h) vanishes 
on [— 1 + cx, 1 — ci] outside $00, a set of measure of order A. 

Finally, arguing as in the end of section 4.5.3 in [3], one proves that the number h satisfies the 
required Diophantine condition if e belongs to the set D described in Theorem 1.6. 
4- Coming back to E. The above results imply that, for e G D sufficiently small, the Lyapunov 
exponent for the equation family (0.1) vanishes on the interval E([— 1 + c\, 1 — Ci]) outside a set of 



the measure m 





dE 


1 


dip 



dip. By the Cauchy estimates for analytic functions and the first step, on 



[-1+ci, 1-ci], we get f- = 0(M(A)). So, m = o(M(A)) = o(|/„|) (which is obtained using (4.56)). 
As Ci can be taken arbitrary small, this completes the proof of Theorem 1.6. □ 
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5. Possible spectral scenarii for small t 



Consider p defined by (1.31). Being expressed in terms of the tunneling coefficients, p is typically 
either exponentially large or exponentially small. In both cases, the results on the spectral properties 
of (0.1) for small r take a very explicit form. In this section, we assume that, for some positive 5 P , 
one has either 

(5.1) p<e- 5 "/ £ , i.e., S h (E)<mm{S v , v (E)}-5 p , 
or 

(5.2) p>e 5 '/ £ , i.e., S h (E) > min{ S V U (E)} + 8 p: 

and describe the spectral results for these two cases. These results easily follow from Theorems 1.3- 
1.6. So, we only outline the analysis omitting most of the elementary calculations. 



5.1. Preliminaries. We assume that e is sufficiently small so as to be in the case of Corollary 4.3. 
The functions £o and £„■ being linear in E, one can explicitly describe the intervals defined by (1.26). 
For sake of definiteness, we assume (4.28). We denote then the leftmost interval defined in Proposi- 
tion 1.1 by /„- and the rightmost by Iq. We use the notations 

I, = and I =[Efr,E%*\. 

Define also 

2 VlTol \1k\J r \l0ln\ 

7o and 7^ being defined in (2.18), and e 1— > f(e) is the function introduced in (1.26). Note that, 

~^~tv x £ Tv,o, — t v x eT V)7T . 
Furthermore, as / tends to zero when e — > 0, for sufficiently small e, one has 
(5.4) t± < Cemax u T V)U and t± x e 2 f h . 

Then, we get 

Lemma 5.1. Let > 2(£+ + t~)A and > 2(t+ — t~)A. Then, one has 



ET = E-tt- Jittf - 2At~ + + t+ E™ = E + f+ - J (t+) 2 + 2At~ + A 2 + t 



Efr = E- t+ + ^JittY - 2At- + A 2 + t", Eg* = E + tt + ^(t+) 2 + 2At~ + A 2 + t+ 
Ift~ <2{tt -t~)A ; then 



t 



K = E- t~ - ^{t-y - 2A ■ t+ + A 2 + t 7l ; 

ift~<2(tt + t-)A, then 



E™ = E-t~ + Jit-f - 2A ■ t+ + A 2 + t~. 



This lemma easily follows from (1.26) and (2.16). 
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5.2. Small p. Now, let us discuss the intervals Iq and I n and the spectrum in them for small p. Let 

(5.5) M(A) = mm{e \f¥ h , e 2 T h /A}, D(A) = max{A, e Vf\ t }. 
Level repulsion. In view of (5.4), Lemma 5.1 implies 

Corollary 5.1. For sufficiently small e, under the condition (5.1), one has 

(5.6) |/ | x M(A), m x M(A), and dist (7 ,4) ~ £>(A); 

i/ie interval Iq is to the right of Eq, the interval I n is to the left of E n , and for u, \i E {0, n}, one has 



(5.7) dist 



M(A), if // = v, 
D(A), if /i + v. 



This corollary in particular shows that there is always a repulsion between the intervals Iq and I n . 
The nature of the spectrum. Corollary 1.1 and Theorem 1.6 show that the nature of the spectrum 
depends on the value of the quantity 

(5.8) A = elog(r 2 |eo|+r 2 |e 7r |). 

Corollary 5.1 and the definitions of £o and ^ imply 

Corollary 5.2. Pick v E {0,7r}. Fix c\ > 0. Define A v = Th/T VyU . One has: 

• There exists C2 > such that, for sufficiently small e, if A < e~ Clj/£ A u , then A < — C2 for all 

EeI v . 

• For sufficiently small e, if A > e Cl ^ £ A u , then A = slog (A/ ' A v ) + o(l) uniformly in E E I v . 

Remark 5.1. Let A* = E\[T} 1 . By Corollary 5.2, for A < A*, on the intervals Iq and /„., one has 
A < — C < 0. Indeed, for v E {0, 7r}, one has A* < epA u , and p < e~ <5p//£ for some positive 5 p . 

5.3. Large p. Now, turn to the case of large p. For the sake of defmiteness, we assume that T v $ < 
T v n . Note that 



T v n T 



2 



(5.9) ^ x - < e -^ +5 ^ £ « 1 



2 



so that t„ x ef V:7T and < t+ + 1~ x T V:0 <C T V:7T . 

We consider two cases when the formulae of Lemma 5.1 take the simplest form. Simple computations 
imply 

Corollary 5.3. Fix < c < 1. Let the conditions (5.2), (4-27) and (4-28) be satisfied. For 
sufficiently small e, one has 

• if A < ct+, then 



|/ 7r | = 2(t+ + A) + ( £ ?;, 7r ) ) |J | = 2(C-A) + o(er t>iW ), dist (/„, Q x e 
Moreover, E is between E™ and E n , E° nt is to the left of E n , and 

Eo — E 1 *, E n — E x e— — , E n — E° nt = 2t+ + o(eT Vy7T ). 

1+ + 



T ' 



if A > c _1 t+, then 



4t+ + o(eT v>7r ), \I \ x e 2 ^- + sT v , , dist (J , Q = 2(A - i+) + o(eT v ,„). 

44 



Moreover, E is to the right of E™, E° ut is to the left of E v , and 

E - E™ = 2(A - t+) + o(ef VtW ), E n - E° ut = 2i+ + o(ef v>7r ), 

(5.10) E -E™=(l-2 A{tt + t « ) ) g where g~e^, g > 0. 

As c in Corollary 5.3 can be fixed arbitrarily close to one, it gives quite a complete picture of what 
happens in the case of large p. The parameter A governing the nature of the spectrum (see (5.8)) is 
now described by 

Lemma 5.2. In the case of Corollary 5.3, one has 

• if A < cty, then 

A|/„ = elog(l + s v p 2 ) + 0(e) where s u = v for v £ {0,7r}; 

I v\ 

• if A > c~H+, then 

Ak = ebg (5^+0(6), and X\ Io = slog ^A + + 0(e) 

6. Periodic Schrodinger operators 

The aim of this section is to give the precise definition of 9 n and to study its properties. Recall that 
A n = \(9 n + 6~ x ) is responsible for the gap between the "interacting" intervals Iq and /„- in the case 
of small p. 

The constant 9 n is defined by the periodic Schrodinger operator (0.2). First, assuming that V is a 1- 
periodic, real valued, L^ oc -function, we recall well known results on this operator (see [8, 2, 13, 15, 22]). 
Then, we define 8 n and analyze its properties. 

6.1. Analytic theory of Bloch solutions. 

6.1.1. Bloch solutions. Let ip be a nontrivial solution of the equation 

d 2 

(6.1) — —il)(x) + V(x)il)(x)=Sil)(x), iGl, 

dx z 

satisfying the relation ip (x + 1) = X ip (x) for all x £ R with A £ C independent of x. Such a solution 
is called a Bloch solution, and the number A is called the Floquet multiplier. Let us discuss properties 
of Bloch solutions (see [8]). 

As in section 1.1, we denote the spectral bands of the periodic Schrodinger equation by [E\, E2], 
[E 3 , E4], . . . , [i?2n+i, E 2 n+2\, ■ ■ ■ ■ Consider <S±, two copies of the complex plane cut along the spectral 
bands. Paste them together to get a Riemann surface with square root branch points. We denote 
this Riemann surface by S. In the sequel, 7r c : S 1— ► C is the canonical projection, and £ denotes a 
point in S and E a point in C. 

For any E £ C, one can construct a Bloch solution i/j(x,S) of equation (6.1). Normalized by the 
condition ip(0,E) = 1, the function S 1— » ip(-,£) is meromorphic on S. All its poles are projected by 
7r c either in the open spectral gaps or at their ends. More precisely, there is exactly one simple pole 
per open gap. The position of the pole is independent of x (see [8]). 

Let : : S 1— > S be the canonical transposition mapping; for any point £ in mathcalS different from 
the branch points, the point £ is the unique solution to the equation n c (£) = E different from £. 
The function x 1— > ip(x,£) is another Bloch solution of (6.1). Except at the edges of the spectrum 

(i.e. the branch points of S), ip(-,£) and ?p(-,£) are linearly independent solutions of (6.1). 
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6.1.2. The Bloch quasi-momentum. Consider the Bloch solution ?p(x, £). The corresponding Floquet 
multiplier £ \— > A (£) is analytic on S. Represent it in the form A(£) = exp(ik(£)). The function 
£ i — > k(£) is the Bloch quasi-momentum. 

The quasi-momentum is defined modulo 2ir. It can be regarded as single valued analytic function 
from S to the cylinder C/(27rZ). Then, the quasi-momenta of ip(x,£) and i/j(x,£) differ only by the 
sign. The imaginary part of the complex momentum vanishes only if ir c (£) belongs to the spectrum 
of the periodic operator. 

Below, we denote the derivative with respect to the local coordinate E = ir c (£) by the dot '. The 
derivative of the quasi-momentum is a single valued analytic function on S. One has 

(6.2) k{£) = -k{£). 

We note that k vanishes only inside intervals projecting onto the spectral gaps. 



6.2. Constants 9 n and A n . These constants are defined in terms of a meromorphic differential on 
S. 

6.2.1. Meromorphic differential Q. On the Riemann surface S, consider the function 

J^ip(x,£) (ip(x,£) — ik(£)xip(x,£)) dx 

(6.3) uj{£) = \ , J - . 

J Q tp(x,£)ip(x,£)dx 

where k is the Bloch quasi-momentum of ip. This function was introduced in [5] (the definition given 
in that paper is equivalent to (6.3)). In section 1.4 of [ J], we have proved: 

(1) the differential Q = uj d£ is meromorphic on S; its poles are located at the points of P U Q, 
where P is the set of the poles of ip(x, £), and Q is the set of points where k(£) = 0; 

(2) all the poles of Q are simple; 

(3) res p n = 1, Wp G P\Q, res q Q = -1/2, Vq G Q \ P, res r fi=l/2, WrePnQ. 



(6.4) Kn (y) = \[6 n (y) + ——\ and 9 n (V) = exp(l n (V)), 



6.2.2. Constants 8 n and A„ . We let 

2^" ( ' y 1 e n (v) 

where 

(6.5) l n (V) = I n(£), 



and g n C S is a simple closed curve on S such that 

• g n is located on C\ a(Ho), the sheet of the Riemann surface S where, for Im7r c (£) > 0, the 
Bloch quasi-momentum of ip(x,£) has positive imaginary part; 

• ^ci^Qn) is positively oriented and going around the n-th spectral gap of the periodic operator 
H . 

The functional l n : V i— » l n (V) is defined on L 2 (T), the vector space of locally square integrable, real 
valued 1-periodic functions. One has 

Theorem 6.1. The functional l n has the following properties: 

(1) it is real valued; 

(2) it vanishes at all even periodic potentials V ; 

(3) for any fixed real s, it is invariant under the transformations V(-) h- > V(- + s). 

(4) /„(■) is non-zero on a dense open set in L 2 (Y); 
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Point 1 in Theorem 6.1 is Lemma 6.1 from [3]. The remaining part of section 6 is devoted to the 

proof of the other points. 

We shall use the representation 

(6.6) l n (V)= l -j Sl(£)-n(£). 



This representation follows from the equality § Q(£) = ~ § g &*(£) which follows from formula (6.5) 
in [3], and the fact that l n {V) £ E. 

6.2.3. The functional l n on even functions. We now assume that V is even and prove that l n (V) = 0. 
In view of (6.6), it suffices to prove that 

(6.7) Sl(£) = Sl{£), 
which follows from the relation 

(6.8) ij)(-x,£) = ij)(x,£). 

Indeed, making in (6.3) the change of the variable x i— > — x and using (6.8) and the relations 
ip(x — 1,£) = e~ lk ^ip{x, £) and ip(x — 1,£) = e lk ^ip(x, £), we come to (6.7). 

Finally, check (6.8). Both functions in (6.8) are Bloch solutions of (6.1). Their quasi-momenta 
coincide modulo 2n. Therefore, these solutions are linearly dependent, and, as ip(0,£) = 1 = ip(0,£), 
they coincide. This completes the proof of the point 2 of Theorem 6.1. 

6.2.4. The translation invariance of l n . Let s be a real number. And let ip(x, £, s) and k(£, s) be the 
Bloch solution and the Bloch quasi-momentum for the potential V s (x) = V(v + s). Then, for £ £ S, 
one has 

(6.9) ^,g, a )= ^ + '' g ' 0) , k(£,s) = k(£,0). 

The proof of (6.9) is similar to that of (6.8), hence, we omit it. Substituting (6.9) into (6.3), one 
sees that 

Q{£,s) = tt{£,0) + d\ogij;(s,£,0), 

where d is the external derivative with respect to £. This, the definition of l n and the first point in 
Theorem 6.1 imply the point 3 in Theorem 6.1. 

6.2.5. For generic V , l n (V) does not vanish. To prove the last point of Theorem 6.1, we use 

Lemma 6.1. One has 

• the functional l n is continuous on L 2 (T) endowed with the natural topology; 

• for any (Vq, V\) £ L 2 (T) x L 2 {T), the function t \— > l n (Vo + tV\) is real analytic. 

Proof. Let Vo £ L 2 (T). Let K be a compact set containing no branch point and no pole of the 
Bloch solution £ \— > ip(x,£, Vq). It suffices to check that 

(1) there exists V, a neighborhood of Vo, such that the functions (£, V) ^ k(£, V), (x,£, V) t— > 
i/)(x, £, V) and (x, £, V) i— > ip(x, £, V) are continuous in (x, £, V) £ [0, 1] x K x V; 

(2) for any £ £ K and x £ [0, 1], the functions t i— > ip(x, £, Vo + tVi), t i— > ip(x, £, Vo + W\) and 
1 1— > k(£, Vq + tVi) are analytic in a neighborhood of 0. 
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This easily follows from the classical construction of the Bloch solutions and the Bloch quasi- 
momentum presented, for example, in [22]. We omit the elementary details. □ 

Now, to prove point 4 of Theorem 6.1, we need only to prove that V i— > l n (V) is not identically 
zero. The facts that l n vanishes on even potentials and that it is translation invariant implies that l n 
vanishes at one gap periodic potentials (as well as on many other "simple" potentials). In section 6.3, 
we prove 

Proposition 6.1. The functional l\ does not vanish identically on the set of two gap potentials. 

This proposition implies that l n does not vanish identically for any integer n. Indeed, let V be a two 
gap potential such that h{V) ^ 0. The potential V n {x) = n 2 V{nx) is also one periodic. Denote the 
corresponding Bloch solution and its Bloch quasi-momentum by if;(x,£,V n ) and k(£,V n ). Clearly, 
one has 

(6.10) ip{x,£,V n ) = i(;{nx,n- 2 £,Vi) and k(£, V n ) = nk(£n~ 2 , Vi). 

The second relation from (6.10) implies that, up to scaling by n~ 2 , the n-th gap for V n (x) coincides 
with the first gap for V] thus, it is open. The definition of to (see (6.3)) and (6.10)) imply also that 
ln{V n ) = h(V)'i hence, l n (V n ) does not vanish. This completes the proof of Theorem 6.1. □ 

6.3. The functional l n on two gap potentials. We now prove Proposition 6.1. Let V be a two gap 
potential. Then, the spectrum of the periodic operator is the disjoint union of two compact intervals 
and a closed half- axis going up to +oo. To the Riemann surface S constructed for this spectrum 
according to the prescriptions in section 6.1.1, we add the point oo. Near to oo, we introduce the 
standard local coordinate r = E" 1 ^ 2 . Then, S becomes compact of genus 2. The idea of the proof 
of Proposition 6.1 is to express the differential Q(£) in terms of standard Abelian meromorphic 
differentials associated to S; for those, the integrals along the closed curve g± can be controlled. The 
proof is naturally divided into several steps. 
The behavior Q at infinity. We prove 

Lemma 6.2. The function uo is holomorphic in a neighborhood of oo and admits the Taylor formula: 

■ 3 rl rx 4 rl o- 5 rl rx 

(6.11) w = / dx dx'{V} + — (V(x)-V(0))dx+ — / dx dx'{V" -V 2 } + o (t 5 ) 

4 Jo Jo 4 Jo 16 Jo Jo 

where, for f(x), a 1-periodic function of x G R, we write {f}(x) = f(x) — f(x) dx. 

We prove Lemma 6.2 in section 6.3.1. It implies that Q is holomorphic near infinity. 
Decomposition of Q in terms of standard Abelian differentials. First, following [16], we recall some 
information on Abelian differentials. The surface S being hyperelliptic of genus two, there are only 
two linearly independent holomorphic differentials on S. Away from the branch points, in terms of 
the local coordinate E = 7T C (£), they have the form 

(6.12) lM£) = -^=, n 1 (E) = -^£— where R(E) = f[(E - E s ). 

Fix Pg5. Let P = tt c (V). Define a meromorphic differential by the formula 

^/R{P) dE 



(6.13) n(£,P) 



E-P' 



The differential £l(E, V) has on S only two simple poles located at the points V and V. For a suitable 
choice of the branch of \/~R, the residues at these points are equal to +1 and —1 respectively. 
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Consider the differential Q(£) — Q(£). As the periodic operator has only two open gaps, the Bloch 
solution ip(x, £) has only two (simple) poles V\ and V 2 located respectively in the closure of the first 
and the second gaps. The description of the poles of Q and the analyticity of Q at 00 imply that 
the differential fi(£) — fl(£) — (Q(£,Vi) + Q(£,V 2 )) * s holomorphic 011 Therefore, there exists 
(Ci,C 2 ) e C 2 such that 

(6.14) n(£) - n(£) = n(£, Vi) + n{E, v 2 ) + c&^s) + c 2 n 2 (£). 

Comparing (6.14) with (6.11), one easily obtains 

(6.15) C 2 = --\ dx dx'{V}(x'), C 1 = - dx dx'{V"-V 2 }(x')--C 2 , a = VE j . 

2 Jo Jo 8 Jo Jo 2 ~l 

Dependence on V\ and V 2 . The two gap potentials with a given spectrum depend only on two 
parameters, namely, V\ and V 2 , the poles of the Bloch solution (see [16]). To prove Proposition 6.1, 
we fix the spectrum and study the dependence of l\ on V\ and V 2 . The representation (6.14) imply 
that 

(6.16) h= <f (n(£,V 1 ,E) + Q(£,P 2 ,E)) + C 1 (V 1 ,V 2 ) <f fii + C 2 (Pi,P 2 ) <f n 2 . 

J 91 Jgi J 91 

Note that the number § Qi is not zero. So, the fact that l\ does not vanish identically will follow 
from 

Lemma 6.3. Let c\ 7^ and c 2 be two constants. The expression c\ d p^p 2 + c 2 d p^ d p 2 does not vanish 
identically. 

Lemma 6.3 is proved in section 6.3.2. This completes the proof of Proposition 6.1. □ 



6.3.1. The proof of Lemma 6.2. First, following [16], we explain how to describe the behavior of ip 
near 00 on S. Consider the logarithmic derivative, x( x ^) = ~ log^(x, £). It is 1-periodic in x\ 
for finite gap potentials, £ 1— > £(•,£) is a meromorphic function on S; in a neighborhood of 00, its 
Laurent series has the form 

(6.17) X = ^-fr-?r 2 + ^r 3 + .... 

The Laurent series for xi x i E) is obtained from that of x{ x i E) by changing the sign of r in (6.17). 
In the neighborhood of 00, the representations of ip and k can be obtained by substituting (6.17) in 
the formulae 

(6.18) v/j(x, £) = exp I i / x{ x i £) dx ) , k{£) = I xi x i^)dx. 



We now derive (6.11). Compute the numerator and the denominator in (6.3). One has 

(6.19) / iP{x,£)<ijj{x,£)dx= / e ^o(x^^)+x^',£))dx' dx=1+ / (v(x) -V(0))dx + O(r A ). 
Jo Jo 2 Jo 
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Furthermore, using (6.18), we also get 



ip(x,£) (i){x,8)-ik{8)x^{x,8)^ dx = i j e iIo(x^',£)+x^,£))dx'_ (j { x ( x ' 7 £)} dx'j dx 

T I! dx Jo dx ^ v ^ x ^ ~ T jf ~ v(0)) dx ~ T (l rfx jf ~ y2}(a;,) 

- F(0)) / rfx'{y}(x')^) +0(r 6 ), 



where the operator {■} was defined in Lemma 6.2. Substituting (6.19) and the result of the last 
computation into the definition of u, we obtain (6.11). This completes the proof of Lemma 6.2. □ 

6.3.2. The proof of Lemma 6.3. We prove the lemma by studying the analytic extensions of the 
functions (Vi,V 2 ) i— > C\{V\iV%) and (Pi,^) C 2 (Vi, V%). The proof is divided into "elementary" 
steps. 

1. Following [16] (chapter II, sections 6 and 5), we describe the set of the two gap potentials giving 
rise to the spectrum [E\, E 2 ] U [E 3 , E4] U [E 5 , 00 [. 

Pick P x and P 2 so that E 2 < P x < E 3 and E 4 < P 2 < E 5 . Recall that R(E) is defined in (6.12). Let 
£ 1— > Ai(£) and £ 1— > A 2 (£) be the solutions to the Cauchy problems 

d\i 1 — — r ( <iA 



2 v / -i?(A 1 ), , — ^ = 2 V / - J R(A 2 ), 

(6.20) <( rf£ v v ; ' and < d£ V V 7 

Ai(0) = Pi (A 2 (0)=P 2 . 

The solutions of these equations are smooth real valued periodic functions of £ £ IR. The range of 
Ai is [E 2 , E3], and the one of A 2 is [E4, E5}. The periods of Ai and A 2 are equal respectively to 

E3 dX , _ r Es dX 

and 



E2 vWOT vW 

One has Hi > S 2 . Set 

(6.21) x= Aa 2 (0-Ai(0)^- 

then a two gap potential giving rise to the spectrum [E\, E 2 ] U [E 3 , E4] U [£5, 00) is described by the 
formula 

(6.22) V(a;) = -2(Ai(0+A 2 (0) + cd, 
where c is a constant depending only on {Ej)i<j<5- 

The potential V is periodic if Si and S 2 are rationally dependent, i.e., there exists / and m mutually 
prime such that ZSi = mS 2 . As Si > S 2 , this can happen only for / < m. 

Assume that ZSi = mS 2 . The period of the potential V is equal to J" ~(A 2 (£) — Ai(£)) d£, where 
S = /Si = mS 2 . For the periodic operator with the potential V, the gaps [E 2 ,E 3 ] and [£4, £5] are 
respectively the /-th and the m-gaps. The points P± and P 2 are the projections on the complex plane 
of the poles of the Bloch solution ip(x,£). The sheets of the Riemann surface where these poles are 
are determined by the choice of the branches of the square roots in (6.20). The poles of the Bloch 
solutions are the parameters indexing the family of the periodic potentials with the fixed spectrum. 
Finally, we note that the potential V satisfies the relation 



(6.23) V 2 (x) - V"{x) = -8 (A?(£) + Ai(£)A 2 (£) + A 2 2 (£)) + ci (A x (£) + A 2 (0) + <*. 
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where C\ and c 2 are constants depending only on (£y)i<i<5- 
In the sequel, we assume that 



(6.24) 



(A 2 (0-Ai(0)df 



i.e., that V is 1-periodic two gap potential whose the first two gaps are open. Note that the condi- 
tions (6.24) are satisfied for a suitable choice of the points (Ej)j = i t ..£. 

We complete this discussion with the formulae expressing C\ and C 2 in terms of Ai and A 2 . Formu- 
lae (6.15), (6.22), (6.23) and (6.21) imply that 

C 1 (V u V 2 ) = 3iF(P 1 ,P 2 ) + c 3 G(P 1 ,P 2 ) and C 2 {V U V 2 ) = iG{P u P 2 ) 

where Pj = 7r c (Vj) for j = {1, 2}, C3 is a constant depending only on (Pj)i<j<5, and 

1 

2 
1 
2 



(6.25) 



P(P,P 2 ) 
G(Pi,P 2 ) 



d£(A 2 



d£(A 2 



Ai)(0 / ^'(Af- A?)(0 
Jo 

Ai)(0 /V(A2-A?)(0 
>/ 



^(a 3 



A?)(e), 

A?)(0- 



2. To study Ai and A 2 for complex P x and P 2 , in (6.20), we consider the 
branch of ^—R{u) analytic in C + fixed by the condition ^J—R(u) < 
for u < E\. We introduce the mapping A : C + 1— > C defined by the 
formula 

AW- ^ 



and study its properties. 

The Christoffel-Schwarz formula (see [19]) shows that A conformally 
maps C + onto D, the domain shown on Fig. 9. In the same figure, we 
have shown how A(X) goes along the boundary of D as A moves along 
the real axis from —00 to +00. Note that 

(6.26) ReA(E 1 ) = ReA(E 2 ) = -2/4, ReA(E 3 ) = ReA(E 4 ) = 2/4. 




A(E 2 ) A{E 3) 
Figure 9: The domain D 



Consider the inverse of A on D i.e. the function £ G D 1— > A(£) G C + defined by the relation 
A(A(£)) = £. It takes real values on the boundary of D. Therefore, using the Riemann-Schwartz 
symmetry principle, one can extend it analytically to the strips S 2 — {0 < Im£ < 21m (y4(P 4 ))} and 
Si = {2Im(A(P 3 )) < Im£ < 0}. In result, we obtain a function satisfying the relations 



(6.27) 
(6.28) 



A(£ + H/2) = A(0 for (gS 2 , 
A(£+~) = A(0 for A(£ + S/2) = AP|) for £ G Si. 

Discuss the singularities of A. Let D be the rectangle {— 2/4 < Re£ < 2/4, lmA(E3) < Im£ < 
1m A(E4)} cut along the real axis from — 2/4 to zero. For £ G D, A satisfies the uniform estimates 

2*i/3(3£) -2/3 j < ^ and | A /(^ + 2e - 5m /3(3^-5/3| < ^-l^ 

(3£) -2 / 3 and £ 1— > (3£)~ 5 / 3 are chosen so that they are both 



(6.29) 



|A(0 



Here, the branches of the functions £ 
positive when £ G M+. 

3. For j G {1,2}, the solution £ i— > Aj(£) to (6.20) is given by the formulae 



(6.30) 



A,(0 = A(£ + 0) where i i = A{P j ). 
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From now on, instead of Pj, we shall use the parameter £y. From (6.30), it is clear, that £ i— > Aj(£) 
is analytic in the strip Sj — £j. 

4. Now, we consider the functions (£1,62) ^ -^(£1,62) and (£1,^2) | — > G^i,^)- To complete the 
proof of Lemma 6.3, it suffices to check that, for sufficiently small positive 8, their partial derivatives 
satisfy the estimates 

(6.31) \dlt 2 F(-5,+5) - F r 8 / 3 | < C5~ 2 and \^ 1$a G(-6,+6)\ < C5~ 2 , 

where F Q 7^ and C > are independent of 5. 

Prove the first of these estimates. Representation (6.25) and (6.30) and relations (6.27) and (6.29) 
imply that 



«9| i6 F(-M) = - /V 2 a? + a;a^k 

J 



~ ' [A&)(A?(0 + A?(£ + S/2)) + (AUO + AU£ + 3/2))Al(0] « 
■a/4 

= - /" U'(£ + ^)(A 3 (£ - i£) + A3(-£ - + (A'(£ - i£) - A'(-f - ^))A 3 (£ + i<f)l d£. 
j -s/4 L J 

Then, by means of the estimates (6.29), we get 

dt^F(-6,+6)=F 8- 8 / 3 + O(8- 2 ), 

POO 

F = -2 • 3- 11 / 3 / {2e m '\t + z)~ 5 / 3 (t - z)~ 2 - e~^ 3 (t - z)^/ 3 (t + i)~ 2 ) dt, 



where the branches of the functions t 1— > (t ± i)~ 5 / 3 are fixed by the condition arg((£ ± i)~ 5 ^ 3 ) — > 
as t — > +00. The constant i*o can be computed by means of the Residue theorem. This gives 
Fo = — 5z7r3~ n / 3 2~ 2 / 3 which completes the proof of the first estimate in (6.31). The second one is 
proved similarly. This completes the proof of Lemma 6.3. □ 
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